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Table A List of scientific terms, concepts and principles used in Unit 16 
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Introduced in earlier Units, Introduced in earlier Units, 
or assumed from general or assumed from general Introduced or developed 
knowledge (GK) Unit knowledge (GK) l Unit in this Unit Page 
T 
absolute luminosity 1 potential energy 3 Bethe cycle 26 
acceleration due to gravity 3,4 potential well 14 black hole 48 
alpha-decay 15 pressure 11 Chandrasekhar limit 38 
angular momentum 4 proton 13 cluster simulation 21 
angular speed 2 Schrödinger’s quantum Coulomb barrier 22 
anti-matter 15 theory 14 Crab pulsar 45 
apparent luminosity 1 special relativity 12 Cyg X-1 52 
axis of rotation 25 spectral line 13 degenerate gas 33 
black body spectrum 1 speed of light 9 & 10 deuteron 13 
centripetal acceleration 2 spin (of electron) 14 Doppler broadening 36 
centripetal force 3 temperature 11 Eddington limit 17 
charge 6 thermonuclear fusion 8 Eddington’s conditions for 
colour of a star j| torque 4 centre of the Sun 19 
conservation of angular tunnelling 15 energy levels in white dwarf 32 
momentum 4 wavelength 9 & 10 event horizon 51 
Coulomb force 6 X-rays 9 & 10, 13 | Fermi energy 32 
de Broglie relation p = h/Agp 14 fusion reaction 23 
derivative 2 general relativity 48 
Doppler effect 1 gravitational bending of light 49 
elastic collision 3 gravitational collapse 11 
electromagnetic spectrum 9& 10 gravitational redshift 50 
electrostatic force 6 Hertzsprung—Russell diagram 6 
electrostatic potential energy 6 kinetic pressure 16 
energy levels 13, 14 main sequence star 32 
equilibrium 4 mantle 41 
frequency 9 & 10 mass-lifetime correlation 28 
galaxy 1 mass limit 
gravitational force 4 for neutron star 48 
gravitational potential energy 4 for white dwarf 37 
Hertzsprung—Russell diagram TV 1 mass-—luminosity correlation 28 
internal energy 11 neutrino 13 
ion 7 neutron production 42 
ionization energy 13 neutron star 43 
kinetic energy 3 Pauli pressure 
kinetic theory of gases 11 of electrons 33 
light 9&10 of neutrons 42 
magnetic field 8 plasma 12 
magnetic flux 8 pulsar 43 
main sequence star TV 1 radiation pressure 17 
mass 3 relativistic softening 38 
mass-energy relation E = mc? 12 Sirius A and B 30 
Maxwell-Boltzmann strength limit for neutron star 48 
speed distribution 1] strong nuclear force 22 
moment of inertia 3 supernovae 41 
momentum 3 thermal conductivity 19 
neutron GK, 15 tunnelling 23 
Newton’s law of gravitation 4 white dwarf 30, 
Newton’s laws of motion 3 TV 16 
nucleus 15 
orbital period 2 
parallax 1 
Pauli principle 15 
photon 13 
plasma 8 
polarization 9&10 
positron 15 | 


‘Following an analogy suggested by Herschel, suppose that an intelligent observer, 
who had never before seen a tree, were permitted an hour’s walk in a forest. During 
that space he would not see a single leaf unfold; yet he could find sprouting seeds, 
small saplings, young, full-grown, and decrepit trees, and fallen trunks mouldering 


back into earth, and in that brief hour he might form a correct idea of the life history 
of a tree.’ 


(From Astronomy by Russell, Dugan and Stewart, 1927.) 


Frontispiece A ‘forest’ of stars. This cluster contains several hundred thousand stars, 
which makes it one of the largest of clusters. It appears in the constellation Centaurus, near 
the Southern Cross, and so is never visible from England. It is named Omega Centauri and, 
at a distance of about 16000 light years, it needs a good telescope to resolve the individual 
stars. Its nature was first recognized by Edmund Halley in 1676 during an expedition to 
St. Helena. Presumably all the stars were formed from one cloud of gas, probably about 
15 billion years ago. 


Study guide 


The relationship of the various components of this Unit is shown below. The study 
times in brackets are a notional guide as to how to split up a total of ten hours; they 
reflect the fact that Sections 1 and 5 are easier than Sections 2, 3 and 4. However 
each Section depends strongly on the previous Sections, so I would not advise 
skipping ahead. If you fall behind, something later may have to be sacrificed. 


The aim of the audio-visual sequences is to take a look at an aspect of a Section from 
a slightly different point of view, rather than to be a link in the chain of argument. 
There are one or two snippets of physics that are not discussed in the Text, but these 
will not be assessed. 


Text Associated tape and television 


Section 1 (14 hours) Band 1: how observations provide the basis for consider- 
ations of stability. (15 minutes) 


Section 2 (24 hours) 
Section 3 (2 hours) TV 16: about white dwarfs. 
Section 4 (2 hours) Band 2: discusses the discovery of pulsars. (15 minutes) 


Section 5* (1 hour) Band 3: black holes and speculation on the future of the 
Universe. (10 minutes) 


One of the main reasons for the choice of subject for Unit 16, apart from its intrinsic 
interest, is that it depends on many aspects of physics developed earlier in the Course. 
This Unit should therefore be helpful for revision, and the time you spend on it 
may be expanded by referring back to earlier Units. To make this easier and more 
profitable, the answers to ITQs and SAQs list explicit back references. Remember 
that Table A also contains references to previous Units. There are no purely revision 
questions in the Unit, however, because the specimen examination paper and the 
Summer School questions fill that role. 


Introduction 


This Course has developed at an increasing pace—from the basic concepts of 
mechanics and electromagnetism, to relativity, and to quantum mechanics. The 
subject for this final Unit—stellar evolution—involves all of these topics. This is a 
subject where we go well outside the scope of common experience, and even outside 
laboratory work, to where the forefront of knowledge merges with speculation. 


The piecing together of the physical picture of stellar evolution must, of course, be 
based on an observational picture. The approaches to some of the observational 
problems—those arising from attempts to measure distance—were described in 
Unit 1, especially in TV 1. We take up the story from there, but the main aim of this 
Unit is to show how stars can be almost stable, and yet slowly evolve. 


Section | is introductory; it gathers together some threads and outlines the problems 
to be solved in understanding stellar evolution. Sections 2, 3, 4, and 5 are about four 
different types of star—main sequence stars, white dwarfs, neutron stars and black 
holes—but all these Sections have the same basic structure and order: 


1 basic physics of a type of star; 

2 prediction of the expected nature of that type of star; 

3 comparison of this prediction with observations; and 

4 further predictions and extrapolations based on confirmed features. 

One might think of main sequence stars as ‘live’ stars, in that they have active energy 


sources inside them, and the other objects as ‘dead’ stars. However, it is the ‘dead’ 
stars that carry the imagination furthest, and take us to unknown areas of physics. 


The order in which the ‘dead’ stars are discussed simply follows the order of in- 
creasing mass. White dwarfs have mass less than about 1.4 times the mass of the 
Sun, and neutron stars have mass between about 1.4 and 2.5 times the mass of the 
Sun, yet this difference is sufficient to give them quite different structures. A small 
increase in mass to above 2.5 solar masses is sufficient to carry us to the totally 
different realm of black holes. 


* Note that $354 TV 16, Shades of black, is relevant to Section 5 of this Unit. 
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1.1.2 


From astronomy to astrophysics 


Stages in stellar evolution 


Introduction 


Unit 1 was mostly concerned with those characteristics of stars that reveal, through 
careful astronomical observations, the distances to the stars. Once the distance is 
known, it is possible to calculate absolute luminosity. A pattern emerges when the 
absolute luminosity of stars is plotted against their temperature, as was first demon- 
strated by Hertzsprung and Russell. In TV 1, The message of starlight we saw that 
the pattern gives a hint that younger stars all have a similar internal structure, and 
that this structure changes drastically as the star evolves. Unit 16 takes up this hint 
and, by applying the mechanics, relativity and quantum theory described in this 
Course, shows how the unseen matter inside older stars behaves in totally unfamiliar 
ways. Only an understanding of physics enables us to uncover this behaviour; we 
can have no direct experience of it. 


The Hertzsprung—Russell diagram 


It is probably a long time since you last looked at Unit 1, so first I want to recall 
those aspects that form the starting point for this Unit. I am going to do this on tape, 
as I only want to refresh your memory, not make a formal exposition. The Frames 
referred to on the tape appear on the next few pages. To make it easier to refer back 
to Unit 1 you will find a brief summary of the tape after the Frames, together with a 
list of references. 


Now listen to the first band on Side 2 of Tape 5. 


Hertzsprung—Russell diagram 
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Summary of Band 1, and of relevant material from Unit 1 


LUMINOSITY The absolute luminosity of a star is the amount of light it gives out; 
its apparent luminosity as seen from the Earth depends on its distance (Unit 1, 
Section 2.1). Few of the 90 nearest stars are included among the 100 stars in our 
Galaxy with highest absolute luminosity. Certainly the nearest star, our Sun, is not 
among the most luminous (TV 1, Television Notes). Ancient astronomers had no 
idea that the Sun is simply a typical star. 


COLOUR The colour of a star is closely related to its surface temperature. The 
wavelength of the peak of a black-body spectrum is inversely proportional to its 
temperature (Unit 1, Section 2.5.1). Many stars have a spectrum that roughly 
approximates a black-body spectrum, though it is normally crossed by many 
absorption lines. To the astronomer the nature of the lines in the absorption spec- 
trum provides a second means for measurement of temperature, and a great deal of 
information on the composition of the atmosphere of the star. 


HERTZSPRUNG-RUSSELL DIAGRAM The plot of absolute luminosity against tem- 
perature shows that many stars are on the main sequence (Unit 1, Section 2.5.3). 
These are young gaseous stars. The radius of a star varies enormously in the later 
stages of evolution. Off the main sequence, stars swell by a factor of 100 to become 
red giants and then collapse to become solid white dwarfs no bigger than the Earth 
(TV 1, Television Notes). The number of stars in an area on the Hertzsprung—Russell 
plot is an indication of the time that stars spend in the corresponding stage of 
evolution (see the quotation on the frontispiece). The rate of evolution depends on 
the nature of the energy source. 


ENERGY SOURCES In a gas cloud contracting under gravitational forces, potential 
energy can be converted into heat. The temperature may eventually rise high enough 
to start thermonuclear reactions (TV 1, Television Notes). However, it is only within 
a fairly narrow range of masses that this gives rise to a star—that is, to a quasi-stable, 
luminous body. 


INTERNAL STRUCTURE As there is no way of observing the inside of a star directly 
(any more than we can see the inside of the Earth or Moon), a picture of a star’s 
internal structure must be built up from calculations based on the equations 
describing the behaviour of the stellar material. For instance, the material of which 
main sequence stars are made behaves like an ideal gas, and obeys the equation PV = 
NkT (Unit 11). This equation, along with a lot of additional information, can be 
used to predict various properties of main sequence stars, and these predictions 
must then be compared with observations. 


Stars evolve very slowly as nuclear reactions change the composition of material 
inside them. Most of the time there is a precise balance between gravitational forces, 
which seek to collapse the star, and gas pressure and radiation pressure, which seek 
to blow it apart. The fineness of the balance is indicated by the following quotation, 
which was referred to on the tape: 


“We can imagine a physicist on a cloud-bound planet who has never heard tell of the 
stars, calculating the ratio of radiation pressure to gas pressure for a series of globes 
of gas of various sizes, starting, say, with a globe of mass 1 kg, then 10kg, 100kg 
and so on, so that his nth globe contains 10”kg of gas. The physical constants em- 
ployed in his calculation have all been measured in his laboratory, and no astro- 
nomical data are required. 


Most of the Table would consist mainly of long strings of nines and zeros. Just for the 
particular range of mass about the thirtieth to thirty-second globes the Table becomes 
interesting, and then lapses back into nines and zeros again. Regarded as a tussle 
between gas pressure and radiation pressure the contest is overwhelmingly one-sided 
except between numbers 30-32. 


No. of Radiation Gas No. of Radiation Gas 
globe pressure pressure globe pressure pressure 
29 0.0016 0.9984 33 0.951 0.049 
30 0.106 0.894 34 0.984 0.016 

31 0.570 0.430 35 0.9951 0.0049 
7 32 0.850 0.150 36 0.9984 0.0016 
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We draw aside the veil of cloud beneath which our physicist has been working and 
let him look up at the sky. There he will find a thousand million globes of gas nearly 
all of mass between his thirtieth and thirty-second globes—that is to say, between 
0.5 and 50 times the Sun’s mass, where the serious challenge of radiation pressure to 
compete with gas pressure is beginning.’ 


(Slightly adapted from Sir Arthur Eddington, Internal Constitution of the Stars, 1926.) 


ITQ 1 A particular star has an apparent luminosity that is 10!+ times 
smaller than the Sun’s apparent luminosity. Its distance is3 x 10° times larger 
than the distance between the Earth and the Sun, and the peak wavelength in 
its spectrum is 800 nm. Plot this star on the Hertzsprung—Russell diagram in 
Frame 4 and deduce the type of star it is most likely to be. 


1.1.3 Stellar formation 


Stars are not born with a functioning energy source. It is generally assumed that the 

first stage in the formation of a star is the contraction, under its own gravitational gravitational collapse 
attraction, of a cloud of gas and dust. The contraction is too slow to observe, but 

there are many gas clouds in the galaxy that are thought to be dense enough to pull 

themselves together. 


Figure 1 shows such an object, silhouetted against the background stars. Its diameter 
is several orders of magnitude greater than that of our solar system, and it is massive 
enough to form several stars. During contraction, the gravitational potential energy 
of the cloud decreases*, and to conserve energy, the atoms and molecules of the gas 
gain speed and are accelerated towards the centre of the cloud. As the density of 
the cloud increases, collisions become more frequent, and the inward motion of the 
particles is converted to random motion—that is, internal energy. The result of this 
is that the temperature of the cloud rises, and a young star begins to shine. But how 
long can this last? Can the gravitational potential energy released in contraction 
account for the huge amount of energy that has been radiated by the Sun? You can 
answer this question by attempting the following ITQ. 


Figure 1 Early stages in star formation. There is a dark area, called the Coalsack, covering 
about 4° near the Southern Cross. It is a cloud of dust and gas, that absorbs nearly all the 
light from stars behind it. Its interest from our point of view is that embedded within it are a 
number of higher density ‘globules’, one of which is shown above. It is about 6 arc minutes 
across (about 1/5 of the angular size of the Moon), which at its distance of about 500 light 
years (5 x 10'*m) means a diameter of about 1016 m. This is considered to be a contracting 
cloud that will heat up and split into a small group of stars. 


ITQ 2 Suppose the Sun has been radiating 102° W for 4500 million years 
(as seems pretty clear from the fact that life on Earth has existed nearly that 
long). What is its total energy output per kg during this time? Compare this 
with the gravitational potential energy released per kg, which you may 
assume is roughly GMs/R, where Mg is the mass of the Sun and Rg is its 
radius. (Refer to the data tables in the Appendix for values of M, and Rg.) 


* Remember, the gravitational potential energy of two objects, with masses m, and m3, 


separated by a distance R, is —Gm,m,/R (Unit 4). In a similar way, the gravitational 
potential energy of a spherical cloud is inversely proportional to its radius. 
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1.2 


It should be clear from ITQ 2 that the gravitational energy released in the contraction 
of a gas cloud down to the size of the Sun is completely inadequate to have supplied 
the energy necessary to keep the Sun at about its present temperature over its whole 
lifespan. Moreover, the Sun is not, so far as we know, contracting at the moment, so 
no gravitational energy is currently being released. 


Although there is still no completely satisfactory model of the energy source, there 
can be no doubt that it arises mostly from nuclear interactions between particles 
in the core of the Sun. As you will see in Section 2.2, nuclear physics plays a central 
part in astrophysics nowadays, and to set the scene for this discussion we must take a 
microscopic view of the material inside a star. 


Stellar matter 


The temperature of the surface of the Sun can be estimated merely from its colour 
(Unit 1) and it is about 5800 K. However, you will see in Section 2 that the tempera- 
ture inside the Sun is much higher; it increases to about 107 K at the centre. As was 
pointed out in Unit 8, at such high temperatures material exists in the form of a 
plasma, a mixture of negatively charged electrons and positively charged ions and 
nuclei. 


The reason that matter is in the plasma state at these temperatures is obvious if we 
consider the energy of the particles. As you may recall from Unit 11, the mean kinetic 
energy <£,;,> of the particles in an ideal gas at temperature T is given by 


[Ekin = 3KT, (1) 


where k is Boltzmann’s constant and has the value 1.38 x 10°77 JK‘, and the 
angular brackets mean ‘the average value of .. .”. 


Thus at 107K, 
CE Nx 138 5610= "7K 102K 
=21 x10 '°J 
= 1300eV 


This is very much larger than the energy required to remove the outer electrons from 
atoms (typically about SeV), and so collisions between atoms will result in one or 
more electrons being knocked out, so that ions are formed. This means that the 
material in stars will be in the form of plasma, rather than in the gaseous state. How- 
ever, although the particles in the plasma are charged, the plasma as a whole will be 
neutral. 


On Earth the most common example of plasma is in a flame, and this type of plasma, 
like that at the surface of the Sun, has a very low density. What is remarkable about 
the plasma in the interior of a star like the Sun is its high density. You can estimate 
the average density of the plasma in a particular star by doing the following ITQ. 


ITQ 3 A planet travels in a circular orbit around a star, with an orbital 
period of one year. The diameter of the star subtends 0.0093 radians as seen 
from the planet. 


(a) Use your knowledge of centripetal forces and Newton’s law of gravita- 
tion to deduce that the mass M of the star is given by 


M = v*7/G 


where v is the speed of the planet, r is the distance between the planet and the 
star, and G is the gravitational constant. 


(b) Use the information about the angular size of the star to deduce another 
expression for the mass M of the star, namely 


M = 42 x 10°7r°p 
where p is the mean density of the star. 
(c) Use the information about the orbital period to deduce that 
v = 2.0 x 1077r seconds! 


(d) Use the results of parts (a)-(c) to deduce the mean density of the star. 


plasma 


1.3 


The star referred to in this ITQ is, of course, the Sun, and the planet is the Earth. I 
find it remarkable that to estimate the Sun’s mean density, one need do no more than 
observe (through very dark glass) that the Sun is just eclipsed by a 4p piece at arm’s 
length! It is typical of astrophysics that the maximum use must be made of every 
observation, to extract all the possible implications. 


The mean density of the Sun is greater than that of water (107kgm 3), and so the 
centre of the Sun must be much more dense than water—about 100 times denser, in 
fact. In such a hot dense plasma there will be an enormous number of collisions, 
and, as we shall see in later Sections, the nature of these collisions holds the key to 
understanding stellar evolution. 


So let us consider what sorts of particles will be present in the plasma at the centre 
of the Sun. We may strike molecules off the list, because the bonds between the atoms 
that make up molecules are far too weak to withstand the collisions. What about 
individual atoms? All of the elements found on the Earth are found on the Sun, 
plus a few extra unstable ones, which have long since decayed away on Earth. How- 
ever, the proportions of the different elements are quite different from on Earth, 
with a far greater preponderance of the lighter elements. 


The lightest element, hydrogen, predominates in the Sun, but in an ionized form, 
that is to say as separate protons and electrons. These particles outnumber all others 
put together by an order of magnitude. When the heavier isotope of hydrogen, 
which is known as deuterium, loses its electron, a nucleus containing a proton and a 
neutron bound together remains. This nucleus is given a special name, the deuteron, 
and it also is found in the Sun. Other atoms will generally be found minus one or 
two electrons or, at the highest temperatures, with no electrons at all. Thus, in the 
centres of stars, collisions between bare nuclei are common and are of prime import- 
ance, as you will see in Section 2.2. 


Of course, there are also enormous numbers of free electrons at the centres of stars, 
having a total negative charge equal to the total positive charge of the nuclei. Besides 
their rather passive role of maintaining charge neutrality they have two important 
effects on the star. First, since they have the same energy distribution as the nuclei, 
they add considerably to the pressure. Plasma at a given density and temperature 
has a higher pressure than gas made up from the same particles, but in the form of 
neutral atoms. Secondly, the electrons are very effective at scattering radiation, and 
so they hinder the passage of radiation to the surface. In the white dwarf stars, which 
we discuss in Section 3, a further process intervenes that gives the electron a still 
greater effect. In these stars it is the electrons that govern the star’s structure, with 
the nuclei playing the lesser role of maintaining charge neutrality. 


There are two other particles that we need to consider. The first is simply the posi- 
tively charged counterpart of the electron, usually called the positron, which was 
introduced in Unit 15. Positrons are formed in many nuclear reactions, both natural 
and man-made, but they do not last long, since they are soon annihilated in en- 
counters with electrons. The other important particle is the neutrino, which has no 
charge, and has a mass which is certainly extremely small, and may be exactly zero. 
If the mass is zero, then the neutrino must travel at the speed of light, like the photon. 
However, it has many other characteristics that put it into quite a different category 
from a photon. (For instance the neutrino has a spin which is the same as that of the 
electron, and different from that of a photon.) Neutrinos are formed in many 
nuclear reactions, and they are difficult to detect directly because they interact only 
very weakly with other particles. 


The main properties of these particles are listed in the data tables in the Appendix. 
You should already be able to recall most of the information about the particles 
listed there, but there is no need to remember the precise values of the masses. 


Summary of Section 1 


1 In seeking to understand the structure and evolution of stars, it is necessary to 
start from two ends: observing stars, and trying to formulate a mathematical model 
for stars. The aim, of course, is to see if the two are consistent. 


2 The quantities that can most readily be measured for a star are the luminosity 
and the (optical) spectrum, which is related to the surface temperature of the star. 
These are the parameters used on a Hertzsprung—Russell diagram. 


deuteron 


neutrino 
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2.1 


3 The mathematical model follows the evolution of a cloud of gas (predominantly 
hydrogen and helium), pulled together by gravity. At first, gravitational energy is 
released in the contraction, and this energy is converted into internal energy by 
collisions of the gas particles. The temperature of the cloud rises, producing a hot 
young star. The temperature is high enough to ionize the gas, producing a plasma of 
ions, nuclei and electrons. 


Self-assessment questions for Section 1 


SAQ 1 Imagine a spherical cluster of stars that covers 1° in the sky and whose 
centre is 1000 light years away. What is the ratio of the inverse squares of the distance 
to the nearest and furthest members? What error is made in ignoring the spread of 
distances and assuming that apparent luminosity is proportional to light output? 


SAQ 2 A ‘white-hot’ iron ball is allowed to cool, while measurements are made of 
its luminosity, colour, radius and temperature. Sketch the path of the measurements 
on a graph of radius R versus temperature T. How would the corresponding measure- 
ments look on a Hertzsprung—Russell diagram (i.e. a plot of luminosity versus colour 
of the ball, rather than of a star)? Remember that the convention is to have red (i.e. 
lower temperature) on the right and blue on the left. 


SAQ 3 What measurements would need to be made, and what arguments would be 
used, to show that a particular star was a main sequence star? 


SAQ 4 Consider hydrogen and sodium placed in the flame of a Bunsen burner. 
Comment on whether a magnetic field applied across the flame would, in principle, 
have any effect on the flame. 


Physical processes governing stellar 
evolution 


If we could make a speeded-up movie of our Galaxy, with say a million years com- 
pressed into a second, we would see the Galaxy make one complete turn in about 
four minutes, with the spiral arms dissolving and reforming as they turned, yet 
keeping their general pattern. Some of the brighter stars would soon explode and 
fizzle out, not even lasting one turn. Dimmer, less massive stars—like our Sun— 
would shine steadily for many turns. Indeed the Sun ranks among the steadier 
sources, with a nearly constant energy output over the long periods required for 
biological evolution. Life on Earth was confined to the oceans for ten turns of the 
Galaxy, and it took another five turns to produce large animals. The dinosaurs 
disappeared a quarter of a turn ago—sixty seconds of the movie. Man came into 
prominence from one frame of the movie to the next. All this time, the Sun was 
providing just the right conditions for the finely tuned reactions of biochemistry. 
And yet even the Sun’s steadiness is superficial; more apparent than real. There has 
been evolution in the interior of the Sun, hidden deep below a steady surface. Rather 
like a chrysalis, the Sun makes big internal organizational changes behind an almost 
unchanging exterior! 


How can we start to apply physics to such complicated situations? Should we say 
that stars such as the Sun seem to be so fixed and stable that the forces must be in 
static equilibrium, so that the ideas contained in Unit 4 apply? Clearly this cannot be 
adequate, since the equations of static equilibrium can say nothing about the 
evolution, albeit very slow, of a system such as a star. 


To analyse the evolution, it is necessary to find the equations expressing the rates of 
change of the properties of the stars. As the evolution is slow, we may hope that it 
would still be valid to use results of calculations that assume static balance. For- 
tunately, this is indeed a good approximation for most stars. In cases for which this 
simplification cannot be made, such as exploding stars, the astrophysicist is faced 
with a problem that takes even the biggest and fastest computers to the limit of their 
capabilities. 


Balance in stars 


The first major step in the understanding of stellar equilibrium conditions was made 
by Sir Arthur Eddington (Figure 2). He made this step without knowing that nuclear 
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fusion is the main source of stellar energy and the main factor determining the 


dynamics of stellar evolution. He confined his analysis, therefore, to the statics of 
stellar structures, and was able, on this basis alone, to reach some remarkable 


conclusions about the conditions at the centres of stars. 


Figure 2 Prof. Sir Arthur Stanley Eddington (right) with Einstein. A major figure in 
astronomy for thirty years, his pioneer work on stellar structure led the way to understanding 
that the energy source of stars relies on the conversion of mass to energy, as described by 
Einstein’s relativistic relation E = mc’. Eddington not only produced a formidable amount 
of original work but was also noteworthy in the communication of ideas in physics and 
astronomy, both in popularizing the major discoveries, and in bringing little-known work 
to the attention of colleagues. 


Following this approach, let us first consider the main forces on a (stationary) cubic 
metre of material at a distance r from the centre of the Sun (Figure 3). 


edge of star 


mass M(r) 


Figure 3 A cross-section of a spherically symmetrical star, focusing attention on a cubic 
metre of material at a distance r from the centre. The pressure and density at r may be written 
P(r) and p(r) respectively. The mass within the radius r is denoted by M(r). 
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Gravitational force 


If we mentally exclude the material outside the radius r, we may think of our cubic 
metre of material, of mass m, as floating on the surface of a sphere of radius r. The 
mass of the material inside the sphere obviously depends on the radius r, and we will 
emphasize this fact by writing the mass of the sphere as M(r). The cubic metre of 
material will be attracted by the material inside the sphere. Assuming that the Sun is 
spherically symmetrical, the force of gravitational attraction acting on the cubic 
metre of material will be, as Newton showed, the same as if the whole mass M(r) of 
the sphere were concentrated at its centre, and will have a magnitude 
Ha SEn (r) 

r 
But what about the material outside the radius r, which we have ignored so far? It is a 
remarkable fact that we can indeed ignore it. If the Sun is spherically symmetrical, 
the combined pulls of the various parts of the shell outside the radius r cancel out. 
Figure 4 makes the result look plausible, but it needs the mathematics of integration 
to prove it. (Newton actually delayed publishing a paper on gravity until he had 
developed integral calculus for this very purpose.) 


A cubic metre of material sitting on the surface of the Earth is supported by the 
rigidity of the solid material beneath it, but a star is made of plasma and there is no 
such support. There are, however, several other types of outward force that counter- 
balance the inward pull of gravity in a star. 


Figure 4 Only the matter inside radius r contributes to the inward gravitational force. 
For a spherically symmetrical star, the net gravitational force on the cubic metre of material 
at A produced by the matter outside the radius r is zero, so the inward gravitational force 
arises purely from M(r). The reason is that fer each segment of material on one side of A 
there is another corresponding segment on the opposite side which, whilst further away, is 
also more massive. When the calculations are made, the gravitational forces due to the 
material outside radius r cancel out. 


Kinetic pressure 


The most important force counterbalancing gravitational attraction is the force that 
arises from the kinetic pressure of the plasma—the pressure due to the motion of the 
particles. As was shown in Unit 11, the pressure P of an ideal gas is given by the 
equation 

PV = NkT (2) 


where N is the number of particles contained in a volume V. Rearranging this equa- 
tion, we get 


N 
P=—kT (3) 
V 
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kinetic pressure 


The number of particles per unit volume (N/V) decreases as the distance r from the 
centre of the Sun increases, just as the number of molecules per unit volume in the 
atmosphere above the Earth decreases with increasing height. Also, the temperature 
T decreases as r increases because energy is generated at the centre of the Sun and 
flows towards the surface. Both of these changes contribute to a decrease of the 
pressure with increasing radius, as can be deduced from equation 3. 


The change of pressure with radius means that there will be a net force due to the 
kinetic pressure on our cubic metre of material, as you can see from Figure 5. The 
magnitude of the upward force F, on the bottom face of the cube is larger than the 
magnitude of the downward force F, on the top face, because P; is greater than P}. 
So there is a net upward force due to the change of kinetic pressure with radius, and 
the magnitude of the force depends on the rate of change of (N/V) with radius 
and the rate of change of temperature with radius. 


increasing r 


increasing density 
and temperature 


Figure 5 The forces due to kinetic pressure on the four vertical faces of the cube (also shown 
on Fig. 3) are equal in magnitude and cancel out. The forces F ı and F,, however, are not 
equal in magnitude because of the different kinetic pressures P, and P, above and below the 
cube. The difference in kinetic pressure arises from the differences in temperature and 
density, and results in a net outward force. The difference in temperature also produces a net 
outward force due to radiation pressure. 


Radiation pressure 


A second outward force arises from what is known as radiation pressure. As the 
centre of a star is nearly always the hottest part, there is a net flow of radiation out- 
wards. Sunshine originates in the centre of the Sun! On the way out, the radiation is 
scattered off electrons and nuclei an enormous number of times, and the net effect 
of the recoil from all these collisions is an outward force, which can be thought of as 
being caused by radiation pressure. The density of radiation energy, and with it the 
importance of radiation pressure, rises very rapidly with temperature—as the fourth 
power of temperature, in fact. Indeed, it is radiation pressure that poses an upper 
limit to the surface temperature of a star. Too hot a star will blow off its atmosphere 
through radiation pressure. This limit, which occurs at about 50000-100000 K, is 
called the Eddington limit, and we will have occasion to recall its existence in connec- 
tion with supernovae. In the Sun, however, the radiation pressure is much smaller 
than the kinetic pressure (see the Table on p. 10). 


Equilibrium and balance of forces 


If our cube of material is to remain in equilibrium at radius r, then the total force 
on it must be zero, that is 


gravitational force due to forcedueto —— - 
force kinetic pressure radiation pressure 


(4) 


radiation pressure 


Eddington limit 
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What is more, if the whole star is to remain in equilibrium, these three forces must 
balance out at all radii*. This means that the rate at which the gravitational force 
decreases with increasing radius must be the same as the rate at which the sum of the 
forces due to kinetic pressure and radiation pressure decrease. Clearly this is a very 
severe constraint on possible structures for a star. 


In fact, there is one other complication to consider. If the star is spinning with 
angular speed œ, then the cube of material must have a centripetal acceleration of 
magnitude wr towards the axis of rotation, where r is the distance of the cube from 
the axis of rotation. This means that there must be a net force on the cube—a 
centripetal force—that is towards the axis and of magnitude mw7r. Equation 4 
therefore becomes 


gravitational force due to force dueto _ centripetal 
force kinetic pressure radiation pressure force (5) 
(inwards) (outwards) (outwards) (inwards) 


and this is the equation that must be satisfied at all positions in a rotating star. Thus 
there are very strong restrictions on the possible internal structures of stars, and these 
restrictions can be expressed in terms of equations involving the rate of change of 
temperature, pressure, density and composition with radius, that is in terms of 
differential equations. 


Differential equations are frequently used to express relationships between physical 
quantities. Such equations can be used to describe general features, whilst at the 
same time providing an easy way of filling in specific details for a particular problem. 
We need go no further than Newton’s second law of motion for an example. In one 
dimension we may write this as 


dt m 

which is a differential equation. It says that the rate of change of the x-component of 
the velocity of a body with time, i.e. the acceleration in the x-direction, is equal to the 
force component F, acting on the body divided by its mass m. This equation is not 
sufficient on its own to work out the x-component of velocity of the body at a 
particular time, even if the values of F, and m are known. Knowledge of F, and m 
only tells us the rate of change of v,, and all of the lines on the v, versus t graph in 
Figure 6 have the same rate of change of v,, i.e. the same slope. However, if we know 
v, at any one instant, we can select the appropriate line on the v, versus ¢ graph, 
and can determine v, at any other time. 


(to, Vo) 


* There can be oscillations about the equilibrium position. Stars can oscillate as a whole— 
like a jelly or soap bubble—and as mentioned on Bands 1 and 2 of the tape, a study of the 
oscillations can give a clue to the internal structure of the star. 
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Figure 6 An example of the solution of 
a differential equation. If all we knew 
about a mass m were that it is acted upon 
by a constant force component F,, so 
that F, = m dv,/dt, its velocity 
component versus time graph might be 
any line parallel to those shown above 
(for v, < c). If, however, v, for the body at 
a given time is known (such as v, = vo 

at time t = tọ), then a particular line is 
singled out, and this line specifies the 
velocity component v, at any time t. 


To return to stars, the temperature gradient—the rate of change of temperature with 
radius—at a given position in a star depends on the thermal conductivity* of the 
material there and on the amount of energy flowing through that material. The 
thermal conductivity depends in turn on the composition of the material, its density, 
and its temperature. Thus the temperature gradient dT/dr at a certain point in the 
Sun can be calculated if the composition, density and temperature are known at that 
point. This problem is analogous to the much simpler example of the velocity 
component of a body that was referred to earlier. In that example, we only needed to 
know the value vo at one instant of time in order to determine the value at any other 
time. In the case of a star, if we know the temperature at one radius, we can work 
out the temperature at any other radius in the star. 


For the Sun, the information that is required is readily available: we know that the 
surface temperature is about 5800 K, and the radius of the surface is 6.96 x 108 m. 
With this information, it is possible to use the differential equation for the rate of 
change of temperature with radius to work back from the surface to find the con- 
ditions at all other radii. The equations reveal a set of interconnections between the 
various parameters—radius, temperature, pressure, composition, and so on. They 
are, however, far too complicated for the mathematician to find an explicit algebraic 
solution. Before the advent of the computer there were only two ways to take the 
problem further; either make severe simplifying assumptions, with the obvious 
danger that the results would be unrealistic, or deduce some limits to the possible 
values of pressure, temperature, etc. by taking extreme values for some parameters. 
Eddington tried both approaches successfully, but the greatest single step forward 
came from the second approach. 


Eddington found that at the centre of the Sun: 
pressure > 1015 Nm“? 


temperature > 7 x 10°K 


ITQ 4 Find the number of particles per cubic metre in an ideal gas with 
pressure and temperature at the lower limits given by Eddington. Convert 
this to a density by making a plausible assumption about the nature of the 
particles. 


The reason Eddington could only give a lower limit is that he did not know the 
internal composition of the Sun—there is not enough convection to bring material 
up to the surface where it can be identified by spectroscopic analysis (i.e. by matching 
lines in the absorption spectrum of the Sun with lines produced in the laboratory by 
absorption in known elements). Eddington’s calculations highlighted the need for 
better information on the composition, but at the same time his lower limit for the 
temperature at the centre was so high that it pointed the way forward to a method of 
finding out more about the composition. 


The point is that at such high temperatures, nuclear fusion reactions take place 
between the nuclei in the stellar plasma. As a result of the nuclear reactions, the 
abundances of different nuclei change with time. Since the thermal conductivity 
depends on the composition, the temperature distribution in the star will also slowly 
change with time. For most of the life of the star these changes are slow enough that 
there is no need to discard the initial assumption of mechanical equilibrium on which 
the differential equations were based. However, as the equilibrium adjusts to the 
changes in composition, the physical characteristics (such as thermal conductivity) 
of the material of the star change slowly with time; they follow a further set of 
differential equations. Thus there is an intricate but fascinating situation inside the 
star, in which the very slow changes in composition cause related changes in thermal 
conductivity, and hence in temperature distribution, and hence in the central 
temperature. This brings us back to the starting point, because a change in central 
temperature changes the rates of the nuclear reactions. 


Just as we needed to tie the solution of the equations for the temperature distribution 
to a known temperature at the surface, so we need to know the initial composition 
of a star before nuclear reactions start, if we are to tie down a solution of the equations 
for the slow evolution of the composition. How is this initial composition to be 


* Thermal conductivity is a measure of how easily heat flows through a material. Thus 
metals have high thermal conductivities and wood has a low thermal conductivity. 


thermal conductivity 


Eddington’s conditions 
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2.2 


2.2.1 


2.2.2 


found? There are a few stars that can currently be observed in their final stages of 
condensation, so we shall have some idea of their composition when they start to 
shine in a few million years! For existing stars such as the Sun, we have, of course, 
no direct measurement of their initial composition and we must rely on a more 
general argument. For reasons which are outside the scope of this Unit, it is thought 
that the ‘Big Bang’ at the creation of the Universe left behind a mixture of about 75% 
hydrogen and 25% helium, and these values are usually taken as the starting mixture 
in calculations of the evolution of stars. If we make this assumption, we now have 
all the items we need to piece together the history of a star. As we shall see, it is the 
course of the nuclear reactions that links the equilibrium conditions at each instant 
to the progressive evolution of a star. It is to these reactions that we now turn. 


Energy sources in stars 


Introduction 


Eddington’s calculation showed that the temperature in the core of the Sun (and 
indeed most of the stars) must be about 10’ K if it is made predominantly of hydrogen 
and helium in a form that can be regarded as an ideal gas. 


In this Section we examine the nuclear reactions that are possible in such material 
at such temperatures, and find that they can indeed release very large amounts of 
energy. For example, if two 0.5 kg pellets of deuterium could be completely fused to 
helium, more than 5 x 10'*J of energy would be released. For comparison, the 
gravitational energy released in a cloud of gas collapsing to the size of the Sun is 
on average only about 10'' J kg~ +, a factor of 5000 less. But what collisions will give 
rise to fusion reactions, and how can the power output be calculated? 


Particle collisions in a plasma 


As was explained in Unit 11, the atoms in a gas are in motion, and at room tempera- 
ture and pressure they frequently collide. These collisions impose on the atoms a 
particular distribution of speeds, called the Maxwell—-Boltzmann speed distribution, 
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Figure 7 The Maxwell-Boltzmann distribution for 300K, plotted (a) in terms of speed 
(for H,, as in Unit 11), and (b) in terms of energy, as is more convenient for this Unit. Due 
to the wide range of speed (or energy) involved, it is very common to plot this distribution 
with (c) a logarithmic speed scale, or (d) a logarithmic energy scale. The lesson is simply 
that in comparing two Maxwell—Boltzmann distributions, one must be very careful that they 
are plotted in the same way! When plotted in terms of energy, the distribution for a given 
temperature is the same for any particle, independent of mass. 
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and the distribution of speeds in hydrogen gas at room temperature is shown in 
Figure 7a. Now we can present the information in this Figure in a more convenient 
form for our present purposes if we plot the distribution of kinetic energy, rather 
than the distribution of speed, and this is done in Figure 7b. The shape of the two 
distributions is roughly similar, as you would expect from the relationship between 
kinetic energy and speed, which at low speeds (v < c) is simply 


Eşin © mv? 


(This Newtonian approximation will be adequate in this Section, but not always in 
later Sections.) 


The Maxwell-Boltzmann energy distribution applies also to particles in the interior 
of a main sequence star, even though the matter there is in the form of a plasma of 
nuclei, ions and electrons. The relation between motion and temperature is so 
fundamental that it applies to these particles just as it does to atoms. So although it is 
a long extrapolation from gas in a laboratory experiment, to a plasma at tempera- 
tures that have only existed on Earth in the explosions of ‘hydrogen’ bombs*, it is an 
extrapolation that rests on very well-established laws of physics. F igure 8 compares 
the calculated energy distributions of particles at room temperature, and at the 
temperature in the core of the Sun. The energies of the particles in solar matter have 
much higher values, of course. Notice particularly the positions on the two graphs 
of the arrows labelling the ionization energy from the ground state of hydrogen, 
marked £;. For hydrogen at room temperature, the particles nearly all have energy 
below the ionization energy, and so are electrically neutral. But in the Sun, most 
particles have energy above the ionization energy, so that even if neutral atoms do 
form, they will be smashed apart in collisions. This is, of course, what maintains the 
material in the plasma state. 
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Figure 8 (a) The distribution of energies of particles (in the form of plasma) at2 x 107K, 
which is thought to be the temperature in the core region of the Sun. Note the position of the 
ionization energy of hydrogen. (b) For comparison, the distribution of energies of atoms 
at 300K, as in Figure 7d, but now on the same scale as the solar particles shown in (a). 


* At the time of writing (1981), some success had been achieved in fusion research aimed at 
reaching these temperatures in controlled laboratory conditions (recall TV 8). It is probable 
that the generation of useful amounts of energy by these methods will occur for the first 
time by the end of this century. 
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The importance of nuclear processes in stellar interiors is perhaps best seen by 
imagining a few seconds in the life of a high-energy proton under those conditions. 
For instance, we might choose to consider protons with energies that are greater 
than 20 keV. Figure 8a shows that the number of particles with,such energies is small; 
in fact there is a probability of only about 1075 (i.e. a chance of 1 in 10°) that a given 
particle will have energy greater than 20 keV. 


A proton with an energy of 20 keV travels at about 2 x 10°ms~'. The density 
of the plasma in a stellar core might be about 10° kgm, which means that there 
are about 10°? particles per cubic metre. Collisions will be very frequent. In fact, 
if you take the radius of a proton as about 10° '*m and consider the volume of an 
imaginary cylinder, with this radius, which the path of the proton carves out, you 
will discover that each such particle collides, on average, many thousands of millions 
of times a second. 


What will happen in such encounters? The vast majority will be simple elastic 
collisions, in which the kinetic energy of the particles is redistributed between them 
in a random way*. But with such an enormous collision rate we must certainly take 
account of less probable outcomes. Sooner or later the particle we are following will 
have a head-on collision, at high energy, with another high-energy proton. Let us 
for simplicity suppose both particles have 20 keV at the collision, and remember that 
both are protons. Since both protons have a positive charge, there will be an electro- 
static force of repulsion between them which slows them down as they approach each 
other, converting their kinetic energy into electrostatic potential energy. We have 
selected the case of a head-on collision between identical particles because both 
particles will come to rest instantaneously at their distance of closest approach and 
this simplifies the calculations. 


The distance of closest approach is of the order of 107 '* m (see SAQ 6). This distance 
is of great significance in physics, because it marks the frontiers of a new realm. At 
distances less than about 10~'*m, a new force comes into play, which attracts all 
nucleons together. This is the strong nuclear force, which you met in Unit 15 in the strong nuclear force 
context of the -decay of nuclei. At small enough distances, the strong nuclear force 
overwhelms the Coulomb force. (The attractive gravitational force is far weaker 
than either.) If we take the particular case of two protons, we find that the repulsive 
Coulomb force (Figure 9a) and the attractive strong nuclear force (Figure 9b) 
balance out at a distance of about 107 !*m. Thus at this distance, the net force is 
zero (Figure 9c), and the total potential energy (Figure 9d) has a maximum value. 


This ‘hump’ in the potential energy is called the Coulomb barrier. Coulomb barrier 
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Figure 9 Forces between protons. (a) The repulsive Coulomb electrostatic force. (b) The 
attractive strong nuclear force. (c) The sum of (a) and (b); the Coulomb force dominates at 
large distances, but is overwhelmed by the strong nuclear force at distances less than about 
10~** m. (d) The potential energy corresponding to the force in (c). To calculate the net 
effect of the forces it is essential to use quantum mechanics, which is expressed in terms of 
potential energy rather than force (Units 14 and 15). Graphs (c) and (d) are related by the 
expression F, = —dE,,,,/dr (Unit 6). 


* The redistribution is random because it depends on the distance of closest approach in the 
collision and this can have any value. It is this incessant random resharing of energy among 
billions of particles that leads to the distributions in Figure 8. 
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The potential energy corresponding to ‘the top of the Coulomb barrier for two 
protons is about 1 MeV. If the combined energy of two colliding protons is greater 
than this value, they will be able to approach closely enough for the net force to be 
attractive, holding the two protons together. This opens up the possibility of a 
reaction in which the two protons fuse together, releasing energy and becoming 


confined in the potential well as if they had fallen into it. It is clear from Figure 9d - 


that if two protons fuse together* to form a deuteron, a large amount of energy will 
be released. As you will recall from Unit 14, particles confined in a well can only have 
discrete energies (Figure 10). The ground state energy for the proton in a deuteron 
is near the top of the well, and, in fact, about 1 MeV of energy is released when two 
protons fuse together to form a deuteron. (Because of the release of energy this 
type of nuclear reaction is of great importance. Most of our current energy supplies 
are derived indirectly from solar energy and therefore from fusion reactions, and our 
future energy supplies may well depend directly on controlled fusion reactions.) 


If colliding particles had to have 1 MeV or more in order to get over the barrier 
shown in Figure 9d, fusion reactions would be extremely rare in plasmas at a tem- 
perature of 2 x 10’K, because the probability of particles having energy in this 
range is extremely low. (The curve in Figure 8a falls exponentially with increasing 
energy, and reaches about e °°°, or 10° 77°, by 1 MeV!) However we must not 
picture the effect of this energy barrier on the protons in classical terms. It is not 
adequate to picture a proton as a point mass or tiny hard sphere. At small distances 
the behaviour of the colliding protons can only be satisfactorily analysed using 
quantum mechanics. One aspect of this is that tunnelling can allow a proton to 
appear on the inside of the potential barrier even though, in classical mechanics, it 
does not have enough energy to surmount the barrier. Figure 11 may help you to 
visualize the process. 
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Figure 11 A proton tunnelling through a barrier. The top half of this Figure shows a 
quantum wave, governed by the potential shown underneath. Do not take this diagram too 
literally. It is only intended to indicate that in the collisions between nuclei in stellar cores, 
the proton wavelength is so long, compared with the size of the nucleus, that the process must 


be pictured in quantum terms, not as a collision between hard spheres. 


* In fact, as we shall see later, when they fuse together a positron and a neutrino are given 
off, so that the reaction is p + p > d + e* + v. The simplest reaction in nuclear terms is 
p + d > #He + y but this would have been a more complicated example in terms of the 
dynamics because two particles of unequal mass do not come to rest in a head-on collision. 
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Figure 10 Proton energy level in the 
deuteron potential well. 
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As was described in Unit 15, it was Gamow, and also Condon and Gurney, who first 
realized the importance of tunnelling and used it to explain «-decay. It was im- 
mediately exploited by Cockcroft and Walton to cause the first artificial nuclear 
reactions, using accelerating voltages provided by a machine that, by classical 
mechanics, was not powerful enough for the job. By the same token, tunnelling is of 
great importance in the Sun, because the temperature of 2 x 10’ K at the centre of 
the Sun would not be high enough to produce a significant rate of fusion reactions 
if nuclei behaved according to classical mechanics. Fusion is possible because the 
tunnelling process waives the need for an energy high enough to surmount the 
barrier. Thanks to the tunnelling process, protons with energies as low as 20 keV 
can undergo fusion reactions at a significant rate. Such energies are available in a 
plasma at 10’ K. Proton fusion reactions are, in fact, believed to be one of the most 
important sources of energy in main sequence stars. 


Nuclear fusion reactions can take place between particles at 2 x 10’ K, thanks 
to the tunnelling process. Fusion works in stars, but only just! 


As fusion reactions proceed, creating new nuclei, the composition of the stellar 
material gradually changes, with the proportion of heavy nuclei slowly increasing 
as they are built up by fusion of lighter nuclei. As we have seen, the equilibrium 
conditions depend on the composition, so the radius and surface temperature (and 
therefore luminosity) of the star must slowly change as fusion proceeds. This in turn 
will affect the temperature in the core of the star, thereby affecting the energy 
distribution of the particles and the rate of the nuclear reactions. As long as light 
nuclei are abundant, the changes are very slow, taking thousands of millions of 
years. But the course of the nuclear reactions, and the physical state of the whole star, 
are indissolubly linked together, and when the supply of light nuclei runs out, much 
more rapid changes are inevitable. 


The evolution of a star is governed by the nuclear physics of the matter it 
contains. 


Fusion reactions in stars 


We have so far mentioned only the simplest fusion reactions, and we must now be 
more systematic in sorting out which reactions may be important in the evolution 
of stars. 


It will be helpful to recall first the relation between the two ways of analysing the 
energy of a nucleus. In the last subsection we considered the variation of potential 
energy with separation for two particles in a collision, because this shows most 
clearly the vital importance of tunnelling. It also shows, in a qualitative way, that 
when the two particles are ‘fused’ together in a nucleus containing both of them, the 
excess energy must be released in some way—often in the form of y radiation (that is, 
electromagnetic radiation of very short wavelength). 


As we do not know enough about the strong nuclear force to calculate the energy 
released in a fusion reaction, it is fortunate that there is a second, more practical, 
way of determining the amount of energy released. Most of the stable nuclei have 
had their masses measured with extremely high accuracy using mass spectrometers*. 
This allows the energy release in a fusion reaction to be found by calculating mass 
differences, and using the relation between mass and energy discussed in Unit 12, 
namely 


AE = Amc? (6) 


* This is not quite so good as it seems! There are more unstable nuclei than stable nuclei, 
and some of these live just long enough to be important in stellar fusion but not long enough 
for their mass to be measured in the laboratory. 
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If nuclei of masses m, and m3, fuse together to form a stable nucleus of mass M, then 
M is always somewhat less than (m, + m3). The mass difference is 


Am = (m, + m) — M (7) 


and this corresponds to an energy Amc?. The important point for us is that measure- 
ments of masses allow us to calculate how much energy is released in a fusion reac- 
tion, and, in many cases, this sidesteps any need for a knowledge of the detailed 
structure of the nucleus. 


As an example, we will calculate the energy released in the fusion reaction* 


proton + proton — deuteron + positron + neutrino 
P=: (p= ad eS S 
1836 1836 3669 1 0 


The numbers underneath the reaction are the masses of the particles in units of the 
electron mass rounded to the nearest unit. The change in mass is 


[2 x 1836 — (3669 + 1)]m, = 2m, 
This mass is converted into energy, and the amount of energy released is given by 


AE = Amc? = 2m,c* = 2 x 9 x 107-3! kg x (3 x 108ms_‘)? 
= 1.6 x 10°13J = 10°eV = 1 MeV 


So the fusion of two protons according to the reaction shown above releases about 
1 MeV of energy, and this energy appears as kinetic energy of the deuteron, positron 
and neutrino. It is interesting to compare the energy released in a typical fusion 
reaction with the energy released when an electron is bound to an ion. You can do 
this in the following ITQ: 


ITQ5 Compare the ionization energy of deuterium with the energy released 
when a proton and neutron fuse to form a deuteron (i.e. a deuterium nucleus) 
and a photon. (The ionization energy of hydrogen is 13.6eV.) 


Not all fusion reactions release energy, but it is easy to determine which ones do 
from a list of nuclear masses. (These can be obtained from atomic masses by sub- 
tracting the mass of the electrons. The binding energy of the electrons is usually 
negligible in this context as you have just shown.) One has simply to see which pairs 
of nuclei, when combined together, give a positive value for Am. Such a search shows 
that the energy-releasing reactions are mostly confined to the lighter nuclei. In the 
Universe as a whole, the nuclei of the light elements are more abundant than the 
heavier ones, and hydrogen is the most abundant of all. In young stars, which were 
formed from such matter, there is, therefore, no lack of fuel for fusion reactions, 
and we may correctly expect that in stars the most common reactions are those 
between very light nuclei. 


Because each reaction provides an energy Amc’, which is released when the particles 
fall into the potential well of the nucleus that they form, the overall rate of energy 
output of a given reaction is the product of Amc? with the rate at which collisions 
leading to that particular reaction occur. How can this be fitted into an overall 
picture? To predict the course of evolution of a star, both the changes in composition 
and the energy output resulting from all the fusion reactions must be calculated as a 
function of time, by following the combined effect of all possible reactions. > 


One of the first attacks on this problem was made by Hans Bethe in 1939. The date 
is significant, because immediately prior to World War II nuclear physics was ripe 
for rapid development in the academic and (as it transpired) military fields. 


* To recall the symbols used to represent particles: 
p means proton 
n means neutron 
d means deuteron—the nucleus of a deuterium atom, consisting of a neutron and 
proton bound together by the strong nuclear force 
e* means positron 
e` means electron 
y means photon 
v means neutrino (or antineutrino) 
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After sifting through the data available at the time, Bethe suggested that the principal 
energy source in young main sequence stars is the following set of reactions*: 


p+tpod+e* +v+144MeV 
p + d > 3He + y + 5.49 MeV 


p + 3He > 3He + e* +v + 19.8 MeV 
In the first of these reactions, two protons fuse to form a deuteron. In the second, a 
proton and a deuteron fuse to produce a 3He nucleus, and in the third, a proton and a 
3He nucleus fuse to produce a He nucleus. Together, these three reactions turn 

2 g! 
hydrogen into helium, and they make up the Bethe cycle. Later it was realized that the 
densities are so high that the first reaction can also take place in the form of a three- 
body collision in which two protons and an electron all come together simultaneously, 
that is 
pte +p—>d +v + 2.46MeV 

and that the reaction 

3He + 3He — $He + 2p + 12.86 MeV 
is also important. 


These reactions are likely to be dominant in young stars for three reasons. First, 
the light elements are the most abundant. Secondly, these reactions involve the 
nuclei with the smallest electrical charge; this means that the Coulomb barrier 
separating them is smallest, and so the reactions between them can take place at the 
lowest temperatures. And thirdly, light nuclei move faster and therefore collide more 
frequently. However, the light nuclei will gradually be used up, and the energy output 
from the reactions of the Bethe cycle will then diminish. This means that the finely 
balanced equilibrium will be lost. What happens to the star then? 


The drying up of its fusion energy source will tend to make the temperature of the 
core of the star decrease. But a decrease of temperature leads to a decrease in both 
the kinetic pressure and the radiation pressure that are supporting the star against 
the force of gravitational attraction. So the star contracts, increasing its density and 
releasing gravitational potential energy, which heats the star. The increased tem- 
perature and density result in an increased internal pressure, which prevents a 
catastrophic collapse. But there is no stable equilibrium that can be reached in this 
way since the star is radiating energy at a greater rate than it can produce it. The star, 
therefore, continues to collapse until the temperature of the core becomes high 
enough for a new set of somewhat heavier nuclei to take part in fusion reactions, 
despite their higher Coulomb barrier. Equilibrium is then restored. The star has a 
smaller core radius than before, though the radius of the surface is not necessarily 
smaller. 


There are two chains of reactions involving heavier nuclei that must surely occur, 
though whether others are also important is still a matter of debate. The two generally 
accepted chains are first 


3He + $He > 7Be + y + 1.59 MeV 
{Be +e > {Li + v + 1.37MeV 
3Li + p > 24He + y + 17.35 MeV 

and secondly 
3He + $He > jBe + y + 1.59 MeV 
{Be + p > 8B + y + 0.13 MeV 
8B > §Be + e* +v + 17.46 MeV (half life = 0.77s) 
$Be > 3He + $He + 0.10 MeV (half life = 2 x 107 !°s) 

ITQ6 (a) Show that the net effect of both these cycles as far as the nuclei 


are concerned is to convert 3He + p to $He with the release of about 20 MeV. 
(Treat the reactions as equations.) 


* The superscripts precedingthe chemical symbols are the mass numbers of the nuclei (the 
total number of protons and neutrons) and the subscripts are the atomic numbers (the 
numbers of protons in the nuclei). 
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Bethe cycle 


2.2.4 


(b) Express 20 MeV per 3He nucleus in terms of joules per kilogramme of 
$He and compare this with your results for ITQ 2. 


Of course, the helium nuclei, which feed these reactions, will in turn be slowly 
depleted, so this second phase of a star’s life must also eventually come to an end. 


There are still other reactions, including another set proposed by Bethe, that assume 
great importance if heavier nuclei such as carbon and nitrogen are present and if the 
temperature rises higher. All possible reactions must be considered if the evolution 
of a star is to be understood. That would be far too ambitious for this Unit. However, 
in later Sections of this Unit we will describe what happens when all the nuclear fuel 
is exhausted. 


Predicted properties of main sequence stars 


It has proved very difficult to test the plausible, but nevertheless speculative, ideas 
about which reactions are most important in main sequence stars. Unfortunately, 
there is no way of getting hold of material from the core of the Sun for chemical 
analysis. The only direct link we have with conditions in the core is via the neutrinos, 
which can traverse the outer layers very easily. It is just within the bounds of tech- 
nical feasibility to detect the neutrinos from $B decays. These neutrinos are the most 
energetic ones produced in the reactions described above, and an experiment to 
count them has been in progress in a mine deep below Colorado for a number of 
years. It is now clear that the rate of detection of neutrinos is lower than originally 
expected on the basis of calculated reaction rates. This discrepancy is not yet under- 
stood—there are too many possible explanations! 


Since no other reactions can be detected individually, the only alternative is to 
undertake a huge calculation, in which all known reactions are taken into account, 
and the changes with time of the radii and temperatures of a complete cluster of stars 
are determined. The results of such a calculation can then be compared directly with 
observations of the luminosity and colour of the stars in a cluster whose age is 
known. In other words, the predicted Hertzsprung—Russell diagram for the cluster 
is compared with that drawn up from observation. The results of such a comparison 
are more encouraging than the neutrino experiment. Figure 12 shows the results for 
the cluster M3; it is clear that the calculation reproduces the general trends of the 
observations reasonably well. This agreement is a bridge between the microscopic 
picture of nuclear processes in stars, and the observations of stellar clusters. The 
keystone of this bridge is the Hertzsprung—Russell diagram. 
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Figure 12 The ultimate test of calculations of stellar structures. Theoretical models of 
stellar structure yield predictions of the radial variation of pressure, radius, temperature and 
composition. Unfortunately these cannot be checked separately and in detail because the 
conditions inside stars are not directly observable. However, the predicted behaviour of a 
cluster of stars can be compared directly with the observed Hertzsprung-Russell diagram 
for a real cluster. There are just two problems: it involves an enormously long and compli- 
cated calculation, and there is no precise way of measuring the age of a real cluster! 


Observations of stellar clusters also show that there is a correlation between the mass 
of a main sequence star and its absolute luminosity (Figure 13). The computer 
calculations also predict such a correlation, and this is additional evidence for the 
correctness of the stellar model. The correlation between mass and absolute lumi- 
nosity is strong enough to allow astronomers to use absolute luminosity (which can 
often be measured) to estimate mass (which usually cannot be measured, except for 
binary stars with convenient orbits). 
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Figure 13 The correlation between mass and absolute luminosity for main sequence stars. 
The absolute luminosity of a star depends in a complicated way on many factors. However, 
for main sequence stars the dominant factor is the mass, and this gives an approximate, but in 
practice very useful, relation between mass and luminosity. A rough fit to the data can be 
obtained by taking luminosity approximately proportional to M? for M < Mg, and an even 
higher power of the mass for M > Mg. This strong dependence gives a wide range of lumi- 
nosity for a modest range of mass (note the logarithmic scales on both axes). 


The problem of stellar evolution is like a huge jigsaw puzzle; nuclear astrophysicists 
are at the stage at which they have unambiguously fitted together small groups of 
pieces, but there are gaps in the picture, and there are some spurious pieces on the 
table that do not belong to the picture at all. There are, on the other hand, some basic 
and definite conclusions that can be drawn, and these will lead into the next Section. 
These conclusions are: 


1 The energy supply of an active star must eventually be depleted. 


2 Astar bereft of energy sources cannot support itself against gravitational collapse, 
because a source of energy is required to maintain the kinetic pressure and the radia- 
tion pressure that prevent collapse. 


3 A complicated mixture of nuclei heavier than hydrogen and helium is built up 
during a star’s lifetime. 


If you are willing to accept these ideas, we can end this Section with a more detailed 
conclusion. The results in Figure 13 show that the luminosity increases very rapidly 
with increasing mass. For instance, you can see that a star twice as massive as the 
Sun emits light energy at a rate that is twenty times greater. (Try drawing in the 
horizontal line on Figure 13 corresponding to twice the mass of the Sun.) Since the 
energy output per kilogram is limited (see ITQ 6), it is a fair guess that a more 
massive star uses up its available energy in a shorter time than a less massive one, 
so that a more massive star has a shorter life. Figure 14 shows the results of some 
quantitative calculations of the lifetime on the main sequence of stars with different 
masses: the lifetime clearly decreases very rapidly as the mass increases. 


It turns out that stars with different masses come to the end of their active life (after 
spending some time as a so-called red giant) with different structures and composi- 
tions. There are, in fact, three very different types of dead star, and these will be 
discussed in Sections 3, 4 and 5 respectively. 
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Figure 14 The correlation between the mass of a star and its lifetime on the main sequence. 
As there is no way of directly measuring the time that a star spends on the main sequence, the 
curve shown above is the result of calculation rather than observation, but is firmly based 
on the data shown in Figure 13. Because energy output rises so strongly with mass, more 
massive stars have shorter lives despite their greater reserves of fuel. In fact lifetime is 
proportional to M~?-? approximately. Notice that the lifetime of a star on the main sequence 
is much less than the age of the Universe (~1.8 x 10° years). Thus many stars in existence 
today must either have completed their active lifetime, or have been formed late in the 
history of the Universe, perhaps from matter thrown out of exploding stars. On this long 
timescale, star formation recycles matter in the Universe. It is only on the human timescale 
that stars are models of constancy and immutability. 


Summary of Section 2 


1 Stars on the main sequence maintain a delicate and stable balance between the 
inward force of gravity and the outward forces due to kinetic pressure and radiation 
pressure. 


2 The equations relating the temperature and density gradients to the composition 
of the star were first written down and solved, for particular cases, by Eddington. 
He showed that the temperature in the central core of stars is greater than 107 K. 


3 This temperature is high enough for fusion reactions to take place between 
particles at the high-energy end of the kinetic energy distribution, thanks to the 
quantum—mechanical tunnelling process. The most important reactions are, generally, 
those with the lowest Coulomb barrier, namely those between the lighter nuclei. 
Most of these reactions release energy in fusion, and it is this abundant source of 
energy that gives stars their long life. 


4 Although there is debate about the exact relative rates of all the fusion reactions, 
overall calculations of the course of evolution of stellar clusters reproduce well the 
observed patterns on Hertzsprung-Russell diagrams and the mass—luminosity 
relation. The calculations also predict how the lifetime of a star depends on its mass. 


Self-assessment questions for Section 2 


SAQ 5 _Isthe ratio of the magnitude of the centripetal acceleration to the magnitude 
of the gravitational acceleration greater at the surface of the Earth or at the surface 
of the Sun? (Make the comparison for points on the equators, and use the value of the 
Sun’s rotational period given in the data tables in the Appendix.) 


SAQ6_ Two protons, each having a kinetic energy of 20000 eV (in the laboratory 
frame of reference) approach head-on. 

(a) What is the speed of each one? (Use a Newtonian approximation.) 

(b) Is the Newtonian approximation valid for these protons? 

(c) What will be the potential energy E,,, of the system at the moment of closest 
approach, taking E,,, = 0 at infinite separation? 

(d) If only the electrostatic forces are important, what is the distance of closest 
approach predicted by classical mechanics? 
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3.1 


3.2 


SAQ7 (a) Will the strong nuclear force tend to increase or decrease the distance 
of closest approach calculated in SAQ 6 (albeit by a tiny factor in this case)? 


(b) If the permittivity of free space, ¢9, could be changed, would it need to be 
increased or decreased in order to reduce the probability of fusion in a head-on 
collision between particles of given energy? 


(c) Would tunnelling have much effect on the rate of fusion reactions if the tem- 
perature was much higher, or much lower, than 107 K? Would tunnelling occur in a 
reaction between a neutron and a nucleus? 


SAQ8 Explain, in one or two sentences, the key reasons why (a) the centre of a 
star must be hot, and (b) the temperature of a star depends on the composition 
throughout the star. 


White dwarfs 


Study comment TV 16 is closely related to the material in this Section. 
You will find a synopsis of this programme in the Television Notes. 


Dead stars with mass less than 1.4 times that of the Sun 


In Section 2 we were concerned with stars in the first phases of their lives, when 
active nuclear energy sources in their cores kept them highly luminous. Such stars 
are made predominantly of a plasma of particles, that behaves as an ideal gas. The 
kinetic pressure exerted by this plasma supports the star against gravitational forces. 
We now turn to objects that are different in all respects. As a star’s internal energy 
sources become exhausted, gravity makes it collapse until it can be supported by 
some pressure other than that due to the kinetic pressure of an ideal gas. In this 
Section we will consider what happens to those stars that have mass less than 
1.4Mg, where Mg is the mass of the Sun. In later Sections we will consider the fate 
of more massive stars. 


An absurd star 


In 1834 Friedrich Bessel, mathematician and astronomer, published his discovery 
of a slight waviness in the motion of the bright star Sirius, the Dog Star, in the con- 
stellation Canis Major*. After many years of observation, the period of this motion 
was found to be about 50 years (Figure 15). 


This is not parallax motion, of course, since the period is not one year. However, 


it was easy to explain the new observation by supposing that Sirius was the bright - 


member of a binary pair of stars, that orbited each other with a period of 50 years 
and with the dim member so coo/ that it was too faint to see with Bessel’s telescope. 
Twenty years after Bessel’s discovery the dim member was seen, thanks to the advent 
of photography, and this confirmed the binary star idea. It was given the name of 
Sirius B, and was so dim that it was difficult to measure its temperature. The first 
measurement was not achieved until 1915, and, though extremely inaccurate, it 
showed a staggering result—the surface temperature of Sirius B was in fact higher 
than that of the Sun! ** If Sirius B was kot and faint, then it had to be very small. In 
fact, the average density appeared to be about 10°kgm 3, about 10 tons per cupful. 
Such tiny white-hot stars have come to be called white dwarfs. 


ITQ7 Look back at ITQ 4 and compare your calculation of the density 
of the core of a main sequence star with the average density of white dwarfs. 
Do you think that matter in white dwarfs can be regarded as an ideal gas? 


* The ‘dog days’ are hot days in late July and August, which are supposed to coincide with 
the appearance of Sirius in the dawn sky. 


** Recent measurements, from outside the Earth’s atmosphere, have shown a characteristic 
black-body spectrum that peaks in the ultra-violet region, at a wavelength of about 110nm, 
corresponding to a temperature quoted as 27000 + 6000 K. (Compare this with the surface 
temperature of the Sun, which is about 5800 K.) 
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Figure 15 Sirius and its companion. (a) The position of Sirius in the sky (at the intersection 
of the arrows). (b) The path throughrthe sky as recorded between 1870 and 1940. (c) The path 
reconstructed as an orbit round an unseen partner, with respect to the centre of mass of the 
pair (marked X). (d) The partner (Sirius B) finally revealed. The spikes on the bright image 
(Sirius A) are spurious effects due to diffraction in the telescope. 


Even Eddington, one of the foremost astronomers of his time, simply could not 
believe those results, and left them on one side as an impossible absurdity. However, 
totally new concepts in mathematical physics were emerging in the 1920’s, and some 
physicists were willing to trust them outside the laboratory and apply them to this 
problem. 


Quantum theory applied at ultra-high density 


In 1926 Erwin Schrédinger published his now famous paper, in which he took the 
quantum ideas that had appeared in an ad hoc form in the work of Bohr and de 
Broglie and set them in a clear framework. This theory of Schrédinger’s has been a 
major ingredient of Units 14 and 15. In the same year, R. H. Fowler applied the 
theory to the densely packed electrons in ‘dead’ stars—stars that have used up their 
nuclear energy, and have collapsed to high density. 


Electrons are fairly free to roam through the body of a star, but they are nevertheless 
confined within it, because any current of electrons out of the star would leave a 
positively charged surface which would attract them back by electrostatic forces. 
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Effectively, the electrons are confined in an enormously wide potential well—as wide 
as the star itself. As you saw in Unit 14, the spacing of energy levels in a potential 
well decreases as the width of the well increases, so the energy levels in a potential 
well as wide as a main sequence star must be very closely packed together. 


Indeed, the energy levels in a main sequence star are so closely packed that the re- 
strictions imposed by the Pauli principle are unimportant. This is because the number 
of states whose energies lie within the range of energies covered by the Maxwell- 
Boltzmann distribution appropriate to the temperature of a main sequence star 
(Figure 8a) is very much greater than the number of electrons available to occupy 
those states. In this situation, the probability of any state being occupied is much less 
than one so it does not really matter that the Pauli principle forbids more than one 
electron from occupying the same state. The electrons can then be pictured classi- 
cally. They are distributed over a wide range of energy levels, but have a well-defined 
mean kinetic energy, given by the classical expression derived in Unit 11, 


[Ekin = 3kT (eq. 1) 


The pressure exerted by the electrons can be obtained from the standard result of 
kinetic theory, 


PV = NkT (eq. 2) 


In deriving the results of kinetic theory in Unit 11, the Newtonian expressions for 
momentum (p = mv) and for kinetic energy (E,;, = mv”) were used. We will 
therefore only be able to use equations | and 2 as long as the particle speeds are very 
much smaller than the speed of light (v < c). 


Let us think what happens when an electron gas in a cubical box is compressed, 

under conditions that keep v < c. We can use a result from Unit 14, that the energy 

levels of particles of mass m in a three-dimensional square well of width D are 
given by 

h? 2 2 2 

= 3mp + ny = = n3) 


As the volume V is equal to the cube of the well width D, we can write 
Eyr (8) 


This result also holds for the electron gas in a spherical star, always provided that 
v < c. Thus as a star contracts, the energy levels will get further apart. This means 
that the number of allowed states for electrons within any particular energy range 
gets smaller as the star contracts. You might expect that the electrons would crowd 
into the reduced number of states that lie within the range of energies covered by 
the Maxwell-Boltzmann distribution appropriate to the star (e.g. Figure 8a). 
However, the Pauli principle prevents more than one electron occupying each 
quantum state. So as the star contracts the electrons are forced to occupy states with 
energies very much higher than (3/2)kT (the average energy of a Maxwell—Boltzmann 
distribution), simply because all of the states with energy less than (3/2)kT are com- 
pletely full. And the more the star contracts, the wider the spacing of the energy 
levels becomes, and the higher are the energies of the states that are occupied. 


In a white dwarf, the distribution of the electrons over the allowed energy levels 
is described by the graph in Figure 16. This shows that essentially all of the states 
with energy less than about 10° eV are occupied, and all of the states with energy 
greater than about 10° eV are empty. The energy at which half of the states are full 
and half are empty is known as the Fermi energy (Ep), and in Figure 16 this energy is 
about 10° eV. The width of the transition region that separates the full states from 
the empty states is of the order of 2kT, and for a white dwarf star at 10°K this is 
about 200 eV. Note that the width of this transition region has been exaggerated in 
the Figure—it is actually very much smaller than the Fermi energy (10° eV). Thus 
even if the temperature of the star changed by a factor of ten, there would not be an 
appreciable change in the average energy of the electrons. 


The important point about the distribution shown in Figure 16 is the very high 
energies involved. At a temperature of 10°K, the average kinetic energy of the 
particles of an ideal gas is given by 


(Ekin) = 3kT = 3 x 1.38 x 10773 x 10°J = 107eV 


This is very much smaller than the average energy of particles that are described 
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by the graph in Figure 16, which shows that states with energies up to about 105 eV 
are filled, so the average kinetic energy would be about 0.5 x 105 eV. 


This is the essential point in understanding white dwarfs; electrons in highly com- 
pressed matter have higher energies than they would if they followed a Maxwell- 
Boltzmann distribution corresponding to the temperature of the white dwarf. 
Clearly, therefore, the electrons in these conditions do not behave like an ideal gas. 
Their behaviour is governed by the quantum mechanics of energy levels and the 
Pauli principle, and the electrons are said to constitute a degenerate gas. 
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Figure 16 In the core of a white dwarf star, all energy levels are fully populated up to an 
energy just below the Fermi energy Ep, which is typically 10° eV. ror the purposes of this 
Section the rectangular distribution cutting off at Ep (shown as a dashed line) is an adequate 
approximation. 


Now I am still assuming at this stage that v < c, so it is still correct to use the relation 
(from Unit 11), 


2N 
P= 3 y (Exin> (9) 


As (Ein) is very much greater than (3/2)kT for electrons in a white dwarf, the pres- 
sure is clearly greater than for an ideal gas (which would follow equation 1). It is this 
enhanced pressure in a degenerate electron gas that supports a white dwarf. Because 
the enhanced pressure arises from the Pauli principle, we will call it the Pauli pressure. 
(It is also called degeneracy pressure in some texts.) Note that because the distribution 
of kinetic energies of the electrons in a white dwarf is practically independent of 
temperature (Figure 16), the pressure in a white dwarf is practically independent of 
temperature! This is, of course, quite different from an ideal gas. 


We can now understand what happens to a main sequence star that has mass less 
than 1.4M; after its supply of nuclear fuel is exhausted. Because its core temperature 
then falls, the kinetic pressure (as an ideal gas) is unable to prevent gravitational 
collapse. However, as the star contracts, the electron Pauli pressure eventually 
becomes large enough to balance out the gravitational attraction and a new equi- 
librium is reached. The star is now a white dwarf, with a radius that is about one 
hundred times smaller than the radius of the main sequence star from which it 
evolved, and with a density of about 10®kgm~?. 


The success of the quantum-mechanical model of white dwarfs goes beyond an 
explanation of the high density. Since the temperature does not affect the Pauli 
pressure appreciably, the radius of a white dwarf remains practically constant as it 
cools down. Thus on a graph that shows radius versus temperature, the points 
representing white dwarfs of a particular mass should be spread along a line parallel 
to the T-axis (Figure 17). Such measurements as have been made conform to this 
expectation, but are unfortunately mainly confined to hot white dwarfs, because the 
cool ones are usually too faint to see. Presumably old ‘red’ dwarf and ‘black’ dwarf 
stars are very numerous, but they are too faint to be visible directly and they don’t 
have companion stars whose position they affect to an extent that can be detected, 
at least with the telescopes currently available. 


We will now discuss three other properties of white dwarfs, namely their period of 
rotation, their magnetic field, and the length of time that a star exists in the white 
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3.3.1 


dwarf state. These properties can be predicted, using simple arguments of classical 
physics, from the quantum-mechanical model of a small, high density star. We will 
first consider the rotational period of a white dwarf. 

Further predicted properties of white dwarfs 


When a star collapses to white dwarf size, its moment of inertia, which depends on 
the square of the radius (Unit 3), is greatly reduced. However, the angular mo- 
mentum must retain its original value (i.e. be conserved) if no external torques are 
acting on the star (Unit 4). Let us take this to be the case for the time being, though 
we will have reason to question it later. 


Now, the magnitude L of the angular momentum is the product of the moment of 
inertia J and the angular speed w, 


E = To. 


So if J decreases in the collapse, whilst L remains the same, the angular speed must 
increase, and the period T of rotation must be correspondingly reduced. 


The equation for conservation of angular momentum can be written as 


initial angular momentum = final angular momentum 
(before collapse) (after collapse) 


Ta; = hor, 
where the subscripts ‘i’ and ‘f’ distinguish the initial and final moment of inertia 


and angular momentum. Suppose, for the sake of simplicity, that stars are spheres 
that have uniform density. For a uniform sphere, 


I = 2MR?, 
where M is the mass of the sphere and R is its radius (Unit 3). Thus 


2 2 =: 2 
=MR;0; = MR 0r, 


-a 
wi = R (0) 


Alternatively, since the period T of rotation is given by T = 2z/w, we can write 


which reduces to 


final period _ Te _ (B) aD 


R. 


initial period Si i 
As white dwarfs are about a hundred times smaller than main sequence stars such 
as the Sun, we might expect periods that are about 104 times less than the period of 
the Sun. The Sun’s rotational period is about 25 days and this means periods for 
white dwarfs of only a few minutes. More sophisticated calculations, taking the 
density variations into account, predict periods of a few hours for a typical white 
dwarf of solar mass. We will show how the measured rotational periods compare 
with this prediction in the next Section. 


As the star contracts, the magnetic field changes, and it is easy to estimate the change 
using the idea of magnetic flux. As you may recall from Unit 8, the magnetic flux 
® through an area is defined as the product of the magnitude B of the magnetic 
field, the area A, and the cosine of the angle between B and the line normal to the 
area: 


® = BAcos@ 


For our purposes, we will assume that the magnitude B of the magnetic field has 
the same value all over the star, and for the moment, we will also assume that it is 
everywhere perpendicular to the surface of the star, so that cos 8 = 1. The rule that 
applies to materials which, like the plasma in a white dwarf, are good conductors of 
electricity, is simply that ® remains constant during the contraction. This is not a 
rigorous law, unlike the conservation of angular momentum, but a useful approxi- 
mate rule. It is sometimes referred to as the ‘freezing in’ of the field lines, for reasons 
that are apparent from Figure 18. Since the surface area of the star decreases during 
collapse, it is obvious that the field must increase. To within the accuracy of the 
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Figure 18 The Sun and a typical white 
dwarf—surface and core temperatures 
and an indication of the magnetic field 
(not to scale). The pattern of field lines 
is the same in both cases but the 

closer spacing around the white dwarf 
indicates that the field strength is 
greater than round the Sun. 
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assumptions we have made,* we can write: 


initial magnetic flux = final magnetic flux 


(before collapse) (after collapse) 
Thus BA, = B, As 
or B;4nR? = B,4nR? 
B, RN 
and so T= (2) (12) 
B: \R 


Comparing this result with equation 10, you can see that the prediction is that the 
magnetic field increases by the same factor as the angular speed increases (assuming 
that both angular momentum and magnetic flux remain constant). 


Finally, let us consider the lifetime of a star in the white dwarf state. By counting the 
number of white dwarfs that can be observed, it has been estimated that their lifetime 
is of the order of tens of millions of years. If their lifetime were shorter than this, 
then we would expect to observe fewer white dwarfs, and, conversely, if their lifetime 
were longer we would expect to observe more of them. But how can a white dwarf 
maintain a high surface temperature for 10’ years in the absence of internal energy 
sources? 


There are two main reasons, both of which are quite straightforward. First, there is 
a lot of energy stored in a hot white dwarf, thanks to its very high density, and this 
stored energy can only escape relatively slowly, in spite of the high temperature at 
the surface, because the surface area is so small. It is the small surface area that 
accounts for the low absolute luminosity of white dwarfs. (Remember, white 
dwarfs are faint stars.) Secondly, the outer layer of atmosphere of a white dwarf is 
less dense than the core, and it therefore has a relatively low thermal conductivity. 
You can think of it as a blanket that limits the flow of energy to the surface of the 
star. Inside the blanket, the dense core is much hotter than the surface of the star. 
However, the thermal conductivity of the high-density core is very large, and so the 
temperature gradient in the core is relatively small. As a result, the temperature at 
the centre of the core is much less than in a main sequence star. This is consistent 
with our picture of core temperatures that are too low for kinetic pressure to be 
important. The combination of small surface area and atmospheric blanket means 
that it takes several million years for a white dwarf to cool down enough to become 
unobservable. 


On an astronomical timescale this is a very short interval. A star with the same mass 
as the Sun has a lifetime of about 10'° years in the main-sequence state (Figure 
14), and this is 500 times longer than its lifetime as a white dwarf. The combination 
of short lifetime and low luminosity explains why so few dwarf stars are observable. 
But if the model of the dwarf star is correct, then there are huge numbers of cold 
dwarf stars in our Galaxy and elsewhere in the Universe. This cold dwarf state is 
the ultimate fate of our Sun—in about 5 x 10° years time—and of all stars with 
mass less than about 1.4M,. More massive stars do not share this fate, for reasons 
that will be discussed in Section 3.5. Sections 4 and 5 will discuss what eventually 
happens to these stars. 


Observations of white dwarfs 


White dwarfs are difficult to see because they are so faint, and so it is not surprising 
that there is a paucity of measurements of their periods and fields. However, in 
favourable cases, the rotational period and magnetic field can be determined from 
measurements of the spectrum and polarization of their light. Let us consider first 
the effect of rotation of a white dwarf on the spectral lines that it emits. 


In Unit 1 you saw that the wavelength of spectral lines emitted by stars that are 
moving away from, or towards, the Earth are shifted relative to the wavelengths that 
would be measured for the same lines emitted from sources in the laboratory. This 
wavelength shift is known as the Doppler effect. The wavelengths emitted from 


* In fact we do not need to assume that B is constant or cos @ = 1, only that the field shape 
remains the same, as in Figure 18. However, the result that we would obtain with this more 
realistic assumption is the same as we will derive using the less realistic assumptions of 
constant B and cos @ = 1. 
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stars that are travelling away from the Earth appear to be longer—they are red- 
shifted—whereas the wavelengths emitted by stars that are approaching appear to 
be shorter—they are blueshifted. Thus measurement of the wavelengths of lines in 
the spectra of stars allow astronomers to determine the speed at which a star (or 
galaxy) is approaching or receding from the Earth. 


Now consider what happens to the spectral lines if a star, such as a white dwarf, is 
rotating rapidly, as well as moving away from the Earth. A spectral line emitted from 
the centre of the face that is towards us—point A in Figure 19—will have a wave- 
length shift determined by the speed at which the star is receding from Earth. The 
rotation of the star will scarcely affect this wavelength, since the rotational motion of 
point A is transverse to the line of sight. Point B, however, is receding at a greater 
speed than the star as a whole, and so light from this edge will be more redshifted 
than light from A. Conversely, point C is receding at a lower speed than the star as a 
whole, and so light from C will be less redshifted than light from A. 


Figure 19 The Doppler effect for a rotating star. The eye represents an observer on the 
Earth and the small sections of wave depict the wavelengths of the light as they will be 
received on the Earth. The differences between the wavelengths have been exaggerated for 
clarity. 


Now, even with the most powerful telescope, a star appears as a point of light, and 
it is impossible to distinguish the light coming from different regions of it. Thus the 
result of the rotational motion is to broaden the spectral lines, and the amount of 
broadening is proportional to the speed at which a point on the surface moves 
relative to the centre of the star. To give yau an idea of the amount of broadening 
caused by the rotation, the spectral lines from a white dwarf with a radius of 
6000 km and a period of 1 hour would be broadened by an amount 3 x 1075 
times the wavelength. This means that a spectral line with 2 = 500 nm would have 
an (extra) width, due to the Doppler effect, of 


500nm x 3 x 10°> = 1.5 x 10°7nm 


The Doppler shift due to motion of a white dwarf, and the Doppler broadening 
due to rotation, are shown schematically in Figure 20. 


A second method for determining the period, that may be used even when no spectral 
lines can be seen, depends on a variation of the overall polarization of all the light 
as the star rotates. We will not consider this method in detail as the polarization 
arises from physical processes outside the scope of this Course. 


The spectra from some white dwarfs do indeed indicate that their periods are as 
short as suggested by the angular momentum arguments. However, the majority 
of white dwarfs show larger periods, in some cases up to 100 days. We will consider 
the range of the periods after considering the magnetic field results. 


As well as the broadening caused by rotation, spectral lines from white dwarfs can 
also be modified by high magnetic fields. Each spectral line is split into several 
closely spaced lines by a magnetic field: this is the so-called Zeeman effect discussed 
in Unit 14. The magnitude of the splitting increases as the magnetic field increases, 
and so measurements of the splitting of spectral lines can be used to determine the 
magnitudes of the magnetic field in white dwarfs. The polarization of the light in the 
spectral lines can also be used as an indicator of field strengths. 


The polarizations that are observed indicate that white dwarfs have a wide range of 
magnetic fields. Some do indeed show enormous fields of about 10° T, as predicted 
by the magnetic flux conservation argument in the previous Section. But there is a 
spread down to more moderate values of the magnetic field than would be expected 
from the rather small range of radii for white dwarfs. What is more, there is little 
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Figure 20 The shift and broadening of 
spectral lines of stars due to various types 
of motion. The broadening of a spectral 
line due to rotation can be readily 
distinguished from the overall shift due 
to any motion of the star. 
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sign of the correlation between short period and high magnetic field that we would 
expect from our simple model in which magnetic flux and angular momentum are 
both conserved. This must indicate that some other processes are involved. Because 
most white dwarfs have lower magnetic fields and longer periods than predicted by 
the simple model, one can only suppose that there are dissipative processes that 
drastically reduce the magnetic flux and/or angular momentum of white dwarfs. 
In addition, it is possible that stars, perhaps during collapse, have their own internal 
dynamo processes, (similar, maybe, to the one that produces the major part of the 
Earth’s magnetic field) which could strongly modify the magnetic flux. Remember, 
though, that measurements have only been made on a very small and probably 
atypical sample—a handful of white dwarfs from the millions of dwarf stars that 
must exist in our Galaxy. 


The true picture of a white dwarf is not yet clear, but will certainly be complicated. 
However, there is no doubt about the very high density of these stars, and so we will 
now tackle the question of just how dense and massive a white dwarf can be. 


The limit of a white dwarf mass 


If we take it as established that a white dwarf is supported by the Pauli pressure of 
electrons, whose energies can be predicted on the basis of quantum mechanics, we 
can now tackle the question ‘How massive a star can be supported by the Pauli 
pressure of the electrons?’ For instance, could one have a white dwarf of five times 
the mass of the Sun? 


In Section 3.3 we discussed the behaviour of an electron gas for which the energies 
were governed by quantum conditions, but for which we could still assume v < c. 
If we imagine further compression of this gas, the energies of the electrons will rise 
further, so that some electrons will have speeds close to c. Above a density of about 
10° kg m3, it is no longer adequate to use the Newtonian relation between momen- 
tum and energy. Neither can we use results from the Schrédinger equation, because 
that is also non-relativistic. 


As you know from Unit 15, there is a relativistic quantum mechanics, based on the 
Dirac equation. It would therefore be logical at this point to repeat the analysis 
sketched out in Section 3.3, but using the Dirac equation. To do this would un- 
fortunately go beyond the scope of this Course. However, there is a simpler route to 
the result of such an analysis. This makes use of the expression for the de Broglie 
wavelength, 


h 


~ Kes 


P (13) 
which remains true under relativistic conditions, and also makes use of the fact that 
the allowed levels in a potential well are those for which an integer number of half 
wavelengths fit into the well. 


We again start from the idea that the electrons in a white dwarf are contained in a 
potential well as wide as the star itself. In such a huge well, an electron behaves 
practically like a free electron in the long distances between collisions, so it is 
plausible (but not rigorous) to use the de Broglie relation given above. 


For our purposes, we need to find the pressure exerted by electrons of the allowed 
wavelengths. We will do this in two stages, first relating pressure to energy, and then 
relating energy to wavelength. This will allow us to separate out the effect of the 
relativistic relation between energy and the magnitude of the momentum. 


The relation between the pressure exerted by particles and their energy can be 
expressed in various ways, but the way that will be most convenient for us, because 
it is true even for particles with speeds near that of light, comes from the fact that the 
higher the pressure of a gas, the more energy is needed to compress it by a given 
amount. In fact the relation is simply 


> dE 


P= = 
r (14) 


and you can prove that this is true for a particular case by doing the following ITQ. 
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ITQ 8 Consider a piston of area A that is moved by a small distance Ax to 
compress a gas in a cylinder. How much energy is transferred to the gas? 
Show that P = —dE/dV. 


Suppose we imagine, for simplicity, a star compressed to half its original diameter 
(i.e. one-eighth the original volume) so that the width of the potential well is halved. 
Then all of the de Broglie wavelengths will be halved and so all of the momenta will 
be increased by a factor of two (equation 13). What effect does this have on the 
energy? The answer depends on whether the speeds of the electrons are comparable 
to the speed of light. For white dwarfs of low mass (less than that of the Sun), the 
speeds of the electrons are low enough to use the Newtonian approximation Epin = 
p?/2m (Figure 21). Thus in this case if the diameter is halved, so that the momenta are 
doubled, then we conclude that the energy of all the states would be multiplied by 
four in the compression. As shown by equation 14, this big increase in energy on 
compression makes such a star very ‘hard’. This is the non-relativistic case described 
in Section 3.3 using Newtonian mechanics and the Schrédinger equation. 
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Figure 21 Therelation between energy and momentum for an electron. For a body of mass m, 


Exin = \/p?c? + m?c* — mc?. For low momenta this is approximated by the Newtonian 
expression Ekin = p*/2m, giving the curved line near the origin. For high momenta it is 
approximated by E = pc, giving the nearly straight line to the right. 


However, if we now turn to a more massive star (M > Mg) where the electron 
energies are higher, then the relativistic relation between kinetic energy and momen- 
tum magnitude must be used, and this approaches Epin = pc, as v approaches c 
(Figure 21). The relation between momentum magnitude and wavelength is not 
changed, however, so that the energy rises /ess on compression than for low energies. 
In the extreme case, the energy is only multiplied by two rather than four when the 
dimensions of the star are halved. This means it is easier to compress a more massive 
star. 


Relativity softens the star. 


On these general grounds we can foresee that there will be an upper limit to a white 
dwarf mass, beyond which the star is too ‘soft’ to support itself. Chandrasekhar 
(Figure 22) made the first calculation of this particular mass and found 1.2Mg. 
More recent calculations have given values of about 1.4M,. This upper mass limit 
for white dwarfs is known as the Chandrasekhar limit. Stars of lower mass can find a 
stable configuration as a white dwarf, but stars of higher mass cannot do so. If, asa 
thought experiment, you crammed a mass greater than 1.4M; into a sphere of radius 
about 6000 km, it could not remain as a white dwarf. So what would happen? That 
is the subject of Section 4. 
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Figure 22 The work of Subrahmanyan Chandrasekhar threads through the whole 
development of the understanding of stellar structure. He was born in Lahore, and studied 
first at the University of Madras. He worked on stellar evolution at Trinity College, Cam- 
bridge, under Eddington from 1933 to 1937, but his pioneering ideas on black holes went 


~ unrecognized. He moved to Chicago and wrote a series of monographs on stellar structure 


stellar dynamics, energy transfer and stability that remain classics in the field. He established 
the Chandrasekhar limit to the mass of white dwarfs, as part of his pioneering work on 
degenerate matter in stars. 


Summary of Section 3 


1 When a star has exhausted its supply of fusion energy, it can no longer be 
supported by the kinetic pressure of the ideal gas or by the radiation pressure, and it 
must therefore collapse to higher density. 


2 At higher densities (above about 10° kg m3), a new form of support pressure 
arises. As the star collapses, the number of quantum states with energy less than 
about kT decreases. The Pauli principle restricts the occupation of each state to one 
electron, and so an increasing fraction of electrons are forced to occupy states with 
energies that are much greater than AT. The pressure they exert is therefore greater 
than the kinetic pressure of an ideal gas at the same temperature. 


3 Thanks to this ‘Pauli pressure’, stars can exist in stable equilibrium at densities of 
the order of 10° kgm 3. They are known as white dwarfs. The first indication of 
their existence was the discovery that Sirius B (the faint, massive companion of 
Sirius A) is not a cool star, but a very small hot star. 


4 When a star collapses to such densities, the rotational period would be expected 
to decrease and the magnetic field to rise, both changing by large factors if there were 
no dissipation. Some white dwarfs show the extreme periods and fields that are 
predicted, but they are not typical. 


5 White dwarfs cool slowly and after about 10’ years are so dim that they are 
extremely difficult to see. Thus a typical Hertzsprung—Russell diagram shows few 
dwarf stars, even though they are common in our Galaxy. 


6 The upper limit to the mass of a white dwarf is 1.4M,, and this is known as the 
Chandrasekhar limit. The limit is imposed by the fact that the electron energy rises 
with density and eventually reaches the relativistic region, where kinetic energy is 
proportional to the magnitude of the momentum. The compressibility of such 
relativistic degenerate matter is too high (i.e. the matter is too ‘soft’) to support the 
star. 


S271 UNIT 16 


39 


4.1 


Self-assessment questions for Section 3 


SAQ 9 Consider the following hypothetical experiment. Suppose a set of main 
sequence stars were formed at a certain time, with masses evenly distributed between 
Mg and 5Mg. 


(a) Referring to Figure 14, make a rough estimate of the time that would elapse 
before white dwarf stars began to form (ignore all mass loss processes). 


(b) Would the average mass of the main sequence stars after this time be greater or 
less than the average mass at the start of the experiment? 


SAQ 10 (a) Starting from equation 9, express the pressure of an ideal gas in 
terms of an average momentum of magnitude p and then in terms of an average de 
Broglie wavelength Aap of the particles. (Ignore the difference between <a °> and 
AaB) 7-) 

(b) Suppose the pressure in a white dwarf is 107° N m~? and that there are 2 x 10%4 
electrons perm?. (These are typical values for a white dwarf.) If we assumed (wrongly) 
that the electrons in a white dwarf were in the form of an ideal gas, and we assumed 
that only the electrons were responsible for producing the pressure, what would be 
the electron wavelength? 


(c) What is the kinetic energy of electrons that have the de Broglie wavelength 
calculated in part (b)? Is this kinetic energy greater than the mean kinetic energy of 
the particles in an ideal gas at 10° K, which is a typical internal temperature for a 
white dwarf? 


(d) Would you refer to this star as ‘hard’ or ‘soft’? 


SAQ 11 (a) Ifthe Sun collapsed to the size of the Earth, but retained its present 
surface temperature, by what factor would its absolute luminosity be reduced? 


(b) How long would it take for the ‘collapsed’ Sun to cool by 10*K? Assume, as 
would be appropriate for a white dwarf, that no nuclear reactions are taking place, 
and that the specific heat of the ‘collapsed’ Sun is the same as for water, 
4x 10? Jkg 4K". 


(The power output of the Sun is about 3.8 x 107° W at present.) 


SAQ 12 A sphere of radius 7 x 108m and of uniform density 1.4 x 10° kgm~? 
rotates once every 25.4 days. This sphere is filled with a uniform magnetic field that 
has magnitude 1074 T. It contracts to a radius of 6 x 10°m, with no loss of energy, 
and it retains uniform density and field strength. Find (a) the new rotational period, 
and (b) the magnitude of the new magnetic field. 


Neutron stars 


There are plenty of stars that have masses greater than the Chandrasekhar limit, and 
that cannot therefore support themselves by electron Pauli pressure. What happens 
to them when they finally run out of energy? The answer is that they end up as even 
stranger objects than the white dwarfs that we discussed in Section 3. Stars that have 
masses in the range 1.4M; to about 2.5M, evolve into neutron stars, which are the 
topic of this Section, and the more massive stars become black holes, which will be 
discussed in Section 5. 


Supernovae 


The more massive main sequence stars have hotter cores, partly for the simple reason 
that they have a thicker layer of insulation round the core. At the higher temperatures, 
the nuclei have greater kinetic energies. This means that larger nuclei (with greater 
charges) can overcome the greater repulsive Coulomb force between themselves, 
and undergo fusion reactions. Thus a wider range of fusion reactions is possible in 
these more massive stars. In a long sequence of reactions, larger and larger nuclei 
are gradually built up until nuclei with a mass number of 56, corresponding to iron 
or nickel, are produced. But there the sequence of reactions that liberate fusion 
energy stops. The reason for this is that the nucleons in the iron and nickel nuclei are 
more strongly bound than in any other nucleus. To form nuclei with greater mass 
numbers, in which the nucleons are /ess tightly bound, would require the absorption 
of energy. So at this stage the star becomes energy-bankrupt. Fusion can no longer 
supply energy. 
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The reactions that now take place are ones that absorb energy rather than release it. 
Broadly speaking these reactions are of two types. In the first type, a gamma ray is 
absorbed by a nucleus, which then splits up, perhaps via a sequence of reactions, pre- 
dominantly into helium nuclei. Such reactions absorb energy very efficiently. The 
second type of energy-absorbing reaction is simply the production of neutrinos, 
which, because of their extremely weak interaction with matter, carry their energy 
out of the star unhindered. There are several ways of producing neutrinos, and the 
higher the temperature the greater the rate of energy loss by this process. 


The essential point is that these reactions have the capacity to rapidly reduce the 
temperature of the core of the star. This results in a sudden decrease in the kinetic 
pressure, so that the core begins to collapse. Calculations show that several things 
should then happen in quick succession. The gravitational potential energy released 
as the core collapses raises the temperature above the Eddington limit (Section 2.1), 
so that the flash of radiation, though brief, is powerful enough to drive the mantle— 
the lower density outer layer of the star—outwards into space. Since the mantle still 
contains light nuclei, fusion reactions may be triggered in it by the sudden heating. 
However, the expansion of the mantle soon cools it down and stops the reactions. 


ITQ 9 Using the expression P = —dE/dV (or AE = —PAV) that you 
deduced in ITQ 8, explain why the mantle cools down as its expands. 


The power released from gravitational potential energy during the implosion of the 
core, together with the power from the fusion reactions during the explosion of the 
mantle, leads to an increase in the brightness of the star of several orders of magni- 
tude. Dramatic events like this, called supernovae, have only been seen in our Galaxy 
a few times during the period of recorded history. However, they are so bright that 
they can be easily seen in other galaxies outside our own. Many examples have been 
recorded and measured. Each supernova has a very short life—typically flaring up in 
a few days and fading away steadily, with its light output decreasing by a factor of 
two every few months. But the remnants of these explosions can still be detected, 
and a typical example is shown in Figure 23. 


Figure 23 A very suggestive object! This ring of gas and dust looks very much as though 
it has been thrown out of an explosion. There are quite a few objects like this in the sky and 
they are often called supernova remnants—even when there is no direct evidence for their 
origin. This object is called the Veil Nebula and was discovered in 1784 by William Herschel, 
using a telescope with an 183 inch mirror. (Incidentally he did all the mirror grinding himself. 
His son John later took this and three other telescopes to South Africa, and used them to 
make the first systematic survey of the southern sky.) 


The mechanism of the supernova process has been worked out in detail, and accords 
well with observations. Indeed, supernovae provide a testing ground for data on 
some nuclear reactions. What is more, they provide sources of known absolute 
luminosity, and so supernovae can be used to measure the distances to the galaxies 
in which they are observed (Unit 1). 
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4.2 A giant nucleus 


It is pleasing to be able to explain such dramatic events as supernovae, but solving 
this problem raises a more difficult one about the fate of the core. Given that the 
supernova phenomenon is restricted to massive stars, it may well happen that the 
core is more massive than 1.4Mg, so that it cannot settle down as a white dwarf. Its 
energy sources are, however, presumably exhausted once the burst of gravitational 
energy is over, and the output of light (and other radiation) from supernovae dies 
down in a year or two*. What, then, happens to the core? Being too massive to be 
supported by electron Pauli pressure, it must continue to collapse until some new 
process can supply a support mechanism. We know from Section 3.3 that as the core 
collapses the kinetic energy of the electrons increases, at first simply due to the rise in 
temperature, but later, and much more strongly, because the spacing of the energy 
levels increases as the volume decreases. Eventually, at densities of about 10'2kgm “3, 
the electrons have enough kinetic energy for the following reaction to take place: 


electron + proton (+ energy) yields neutron + neutrino neutron production 
SS E De V. (15) 


This reaction differs from the fusion reactions that we have been considering in an 
important respect, namely that there is no electrostatic barrier to hinder the process. 
Indeed, the electron and proton attract each other. The most important distinction 
is that it is an energy-absorbing reaction, because the products have a greater com- 
bined mass than the reacting particles. 


ITQ 10 Use the information in the data tables in the Appendix to calculate 
the mass difference between the products of the reaction in equation 15 and 
the reacting particles. Hence determine the minimum kinetic energy of the 
reacting particles, if the reaction is to occur. 


If we had taken the data to more significant figures we would find that the electron 
must in fact have a kinetic energy greater than 1.3 MeV for the reaction to occur. 
As you can confirm by doing the following ITQ, electrons with this energy are well 
into the relativistic region, i.e. their speed is comparable with the speed of light. 


ITQ 11 The relationship between total energy and kinetic energy in rela- 
tivistic mechanics is 


total energy = kinetic energy + rest mass energy 


2 


= = 
c? 


Use this relationship to find the speed of an electron that has E,;, = 1.3 MeV. 


Ekin IF me? 


So the neutron-producing reactions described by equation 15 only occur when the 
electrons have relativistic speeds. This means that they don’t occur in white dwarfs: 
remember, for the Pauli pressure to be ‘hard’ enough to prevent gravitational collapse 
in a white dwarf, the electrons must not be in the relativistic region. For stars with 
mass greater than the Chandrasekhar limit, the electrons do have relativistic energies. 
Relativity softens these stars so that the Pauli pressure is inadequate to prevent 
collapse, and it is only in such stars that neutrons are produced. 


The production of neutrons makes a new support process possible, because neutrons, 
like electrons, obey the Pauli principle. The neutrons will be much slower moving 
than electrons because they are so much more massive (m, = 1839m,), and also 
because a large fraction of the energy is carried away from the neutron-producing 
reaction by the neutrinos. It was recognized by Landau and Bohr in 1932, im- 
mediately after the neutron was discovered, that a ‘gas’ of relatively slow-moving 
neutrons could be ‘hard’ enough to make a stable star. They speculated that stars 
might exist composed entirely of the new particles, supported by neutron Pauli neutron Pauli pressure 


* Although no supernovae have been seen in our Galaxy since the invention of X-ray 
detectors, there is indirect evidence for the emission of X-rays from supernovae. These X-rays 
would cause chemical changes in the Earth’s atmosphere, particularly in the concentration of 
oxides of nitrogen. Evidence for such changes has been found preserved in the Antarctic ice 
by drilling down to layers deposited long ago, and making chemical analyses of the samples. 
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pressure, rather than electron Pauli pressure. The neutrons would be closely packed 
together, and they predicted that such a star could support itself if it had a density 
of the order of 10!°kgm >. Indeed a neutron star was pictured as a giant nucleus. 
Not long after this, the suggestion was made that supernovae might produce such 
neutron stars. 


The collapse of a supernova core to a neutron star would be so rapid (a matter of 
days rather than years) that one would expect angular momentum and magnetic 
flux to be conserved, leading to short periods of rotation (of the order of seconds) and 
surface field strengths up to 10°T. In comparison, the contraction of a star to a 
white dwarf is probably far more leisurely with correspondingly more opportunity 
for dissipative processes. Table 1 shows a comparison of the observed characteristics 
of white dwarfs with the even more extreme characteristics predicted for neutron 
stars. 


Table 1 A comparison of white dwarfs, neutron stars and the Sun. 


White dwarfs Neutron star 
Characteristics (observed values) (predicted values) Sun 
Typical mean density /kgm~? 108 10!* to 1018 1400 
Typical mass < 14M, >1.4M, Ms 
Typical radius / km 6000 10 6.96 x 10° 
Typical rotational periods 0.02 to 100 days 0.1s 25 days 
Magnetic field strength/T up to 1000 10° toz 
Moment of inertia / kg m? 3x02 1078 6 x 10*¢ 
Surface temperature / K up to =30 000 ~ 10° 5800 


Pulsars as neutron stars 


At the time of Landau and Bohr’s work, no stars were known with the predicted 
characteristics shown in Table 1. The discovery of neutron stars was long awaited, 
and came finally in an unexpected way. 


As you must have observed, stars appear to twinkle—their apparent luminosity 
fluctuates irregularly on a timescale of seconds. This twinkling is not due to changes 
in the absolute luminosity of stars, but is caused by turbulence in the Earth’s atmos- 
phere. ‘Twinkling’ occurs in the radio waves emitted by stars, too, but this is caused 
by the inhomogeneity of the interplanetary medium through which the radio waves 
from the stars have to travel before reaching the Earth. Anthony Hewish, at Cam- 
bridge in 1964, set out to study the properties of this interplanetary medium by 
observing the radio twinkling of bright sources known as quasars*. His group built 
a radio aerial in the form of a four-acre ‘mattress’ of wires strung out a few feet 
above the ground, as shown in Figure 24. It was a very crude aerial by radio 
astronomy standards, but designed specifically for the job of detecting rapid varia- 
tions in radio signals. Accordingly, the output of the receiver was fed to a recorder 
with a much faster response than normally used at that time with radio telescopes. 
One of the students in Hewish’s group, Jocelyn Bell (now Jocelyn Burnell), found 
that at certain times the telescope picked up regular short pulses of radio waves, 
as shown in Figure 25. Indeed they were far too regular to be explained by 
twinkling; could they perhaps have been due to interference from some Earth- 
bound source of radio waves? 


Study comment Now listen to Band 2 of Tape 5, on which Dr. Burnell 
describes her work. You will need to refer to Figures 24 and 25 in this Unit 
while listening to the audio cassette. 


The following are the key points in the discovery: 


1 It was deduced that the pulses had an extra-terrestrial source because they came 
from a fixed direction in space, that is they were observed by the telescope at intervals 
of 1 sidereal day, not one solar day, (i.e. one Earth rotation with respect to the stars 
rather than with respect to the Sun.) The source of the pulses was called a pulsar. 


2 The pulses and the period were so short that it was concluded that the source 


must be small. 


* Quasars are galaxies with a large power output, in some cases extending to the radio 
region. All known examples are in very distant parts of the Universe. 


neutron star 


pulsar 
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Figure 24 The telescope used in the discovery of pulsars. 
(a) 2048 dipole antennae, making up a four-acre array, were constructed (by graduate 
students) in a field near Cambridge. The telescope was designed for a wavelength à = 3.68 m 
(frequency 81.5 MHz). To give some idea of the effort involved here is a list of the main 
items: 

1.67 tons of copper wire 

14 miles of support catenary 

8.5 miles of coaxial cable 

70 miles of wire for reflecting screens. 
(b) A plan view of the dimensions of the telescope. 
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3 The repetition rate was very accurately constant, which was suggestive of a 
rotating or oscillating source. 


It is now thought that pulsars are rotating neutron stars. The picture would not be 
complete, however, without the link between pulsars and supernovae. This came 
when a pulsar was discovered in the central region of the Crab nebula, which is the 
remains of the expanding mantle of the supernova explosion observed in 
A.D. 1054. This pulsar is exceptional in being detectable by its visible radiation as 
well as by its radio radiation, a fact that was proved in the novel way shown in 
Figure 26. 


Crab pulsar 


Figure 26 The Crab pulsar can be identified by the pulsing of its optical radiation. A 
mechanical shutter in front of a camera was set to rotate at precisely the frequency given by 
the radio observations. The four pictures correspond to different phases of the shutter. 
In (a) the phase is set to block out the light pulses, in (b) it is set to transmit them, and (c) and 
(d) correspond to intermediate phases. The result is a striking identification of the (optical) 
star that is giving the radio pulses. A normal photograph (with a long exposure) would not, 
of course, reveal the variations of the pulsar. 


Further features of neutron stars 


The period of a pulsar can be measured so accurately that there are now many cases 
in which a very gradual increase in period can be detected. The increase in rotational 
period indicates a loss of rotational energy from the pulsar. This energy loss can be 
accounted for in terms of transfers of mechanical and electromagnetic energy 
between a neutron star and the surrounding nebula. A model using a white dwarf, 
with far greater moment of inertia, would not account for the changing period. This 
is the strongest evidence for the identification of pulsars as neutron stars. 


Crab pulsar 
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It is not known exactly how the rotational energy of a neutron star is converted to 
radio waves and other electromagnetic energy, but the observed pulsing of radio 
signals is thought to indicate that radiation is concentrated in a rotating beam 
rather than being isotropic. Thus pulses are received on the Earth as the beam sweeps 
across us rather like the beam of light from a lighthouse. Pulsars do not pulse! This 
beaming could be accounted for if (as with the Earth) the axis of the magnetic field 
were at an angle to the axis of rotation. Imagine the tiny pulsar rotating fast inside a 
large region—known as a magnetosphere—that contains a low density of charged 
particles. The charged particles in the magnetosphere of the pulsar would be sub- 
jected to a varying magnetic field. This creates an electric field (Unit 8), so that the 
particles would be accelerated, and the accelerated charges would produce electro- 
magnetic radiation. Because the magnetic field of a pulsar is so strong, this is 
expected to be a powerful and efficient process. However, the exact details of the 
process are by no means clear. 


The interior of a neutron star may have a complicated structure, since the equilibrium 
proportions of neutrons, protons and nuclei depend on density. Near the surface, 
where the density may be about 10'*kgm_°, there will be a preponderance of heavy 
nuclei. Deeper down, the nuclei will break down into neutrons and some protons, as 
suggested in Figure 27. In the core of the star the density may be about 1018 kgm “3, 
which is even higher than the density of a heavy nucleus. Such conditions cannot be 
created in the laboratory, even during the brief instant of a nuclear collision. 
Theoretical physicists suggest that a ‘gas’ of unstable excited nucleons* may exist 
there. This time no one will say ‘but that’s absurd’! 


excited nucleons ? 


superfluid neutrons 


neutrons 


neutrons and nuclei 


—— magnetosphere 


Figure 27 The internal structure of a neutron star. A neutron star is not all pure neutrons. 
Calculations suggest that there will be nuclei near the surface, and excited nucleons near 
the core. Between them, the bulk of the star will be made of neutrons, but in different physical 
states, according to the density and pressure. The diagram above is based purely on calcula- 
tions, using the results of laboratory experiments on nuclei at high energies. 


Pulsars are a beautiful example of the links between astrophysics and nuclear 
physics. They show how a knowledge of the laws of physics, which we believe to be 
universal, provides a map to guide us to, and through, totally unfamiliar territory. 
The astronomical observations provide, every so often, a signpost to tell us whether 
we are still on the right path. I hope you have confidence in the map and signposts. 
In Section 5, the final Section of both this Unit and the Course, we shall follow them 
as far as anyone has yet been. 


* The term nucleon is used to cover protons and neutrons and also heavier particles with 
similar properties. 
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Summary of Section 4 


1 The collapse of a stellar core, for instance in a supernova implosion, can raise 
the density to a level at which the electrons, through the effect of the Pauli principle, 
have sufficient energy to combine with protons and so produce neutrons. 


2 The neutron energies are low enough, owing to their high mass, to provide the 
‘hardness’ that is characteristic of a non-relativistic gas, so that a stable neutron 
star can be formed. The density corresponding to this state of matter is at least 
10'+kgm 3, which is comparable to the density of an atomic nucleus. One may even 
think of a neutron star as a huge nucleus. The density at the centre of a neutron star 
may be as high as 10'S kgm °. 


3 In 1967, small sources were discovered that emit radio radiation, which is 
received in the form of short pulses. They were soon identified with neutron stars. 


4 The periods of these pulses are of the order of seconds and, assuming these to be 
the rotational periods of neutron stars, are consistent with negligible loss of ro- 
tational energy during the collapse. It seems very plausible that high magnetic fields 
are produced in the collapse, and that these fields generate electromagnetic radiation 
through their action on the magnetosphere of the star. The radiation is generally 
thought to be emitted as a rotating beam, like a lighthouse. This explains why the 
radiation is received in the form of pulses, though the precise mechanism for 
generation of the radiation is still not clear. 


Self-assessment questions for Section 4 


SAQ 13 Why could a star with the parameters given in the middle column of 
Table 1 not be supported by electron Pauli pressure? 


SAQ 14 The discovery of white dwarfs and the discovery of neutron stars both 
led to a long process of scientific investigation, and in both cases careful arguments 
were needed to convince sceptical scientists of the nature of these sources. 


(a) What explanations were put forward initially to account for Bessel’s discovery 
of the path of Sirius, and Bell’s discovery of regular pulses? 


(b) What evidence led to widespread acceptance of current models of these sources? 


(c) How was the development of the theoretical models related to the experimental 
discoveries? 


Beyond nuclear matter 


Does the story of the previous two Sections stand further iteration? Let’s first recall 
the steps in the story. White dwarfs are supported by Pauli pressure of the electrons. 
There is an upper limit to the mass of a body that can be stable in white dwarf 
conditions, given by the effect that we called ‘relativistic softening’—the inadequacy 
of the Pauli pressure when the electrons reach speeds comparable with the speed of 
light. When this upper limit of mass is exceeded, collapse continues past the density 
of a white dwarf until the picture is changed by nuclear reactions that produce vast 
numbers of neutrons. Now try and rewrite the last two sentences changing ‘white 
dwarf’ to ‘neutron star’ and ‘electron’ to ‘neutron’. Does it work? Presumably there 
is an upper limit to neutron star mass that is reached when the neutrons approach 
the speed of light. But when this mass is exceeded, what changes the picture? The 
density of a nucleon is about 10'*kgm~ >, so how can there be nuclear reactions 
leading to a markedly more condensed arrangement? It looks difficult to write the 
second remake of the white dwarf story! 


Modern nuclear physics confirms this feeling that neutron stars are the end of the line, 
and it gives a quantitative limit to the maximum mass of a stable star by posing an 
upper limit to the rigidity of any form of matter. In fact there are two ways to set 
such an upper limit and both are instructive. First, one may consider the way 
Chandrasekhar found the upper limit to the mass of a star in the white dwarf con- 
dition, and see if that method of calculation can be used again for neutron stars. 
Chandrasekhar could calculate a limit because he knew how electrons interact with 
each other via the Coulomb force, and so he could make exact (relativistic) calcula- 
tions of the electron energy levels in a white dwarf and the way they varied with the 
size of the star. It was a good assumption to consider only the electrons, and neglect 
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the nuclei. To repeat this method for a neutron star one needs to know how neutrons 
interact with each other via the strong nuclear force. 


Although much has been learnt about the strong nuclear force from high energy 
particle physics experiments, its effects cannot be calculated as well as for the 
Coulomb force. It does not follow a simple explicit law as does the Coulomb force. 
Moreover it is not only neutrons that have to be considered. There are many other 
types of nucleon that may well be important at the highest densities, and their collec- 
tive behaviour in a neutron star may be different from their individual behaviour as 
studied in collisions. Because of these difficulties, the upper limit to the mass of a 
neutron star cannot be pinned down precisely, but is thought to lie in the range of two 
to three times the mass of the Sun—we will fix on a figure of 2.5 for simplicity. 


The second approach to the problem confirms that there must indeed be an upper 
limit, even though its exact value may be uncertain. It involves some interesting 
and very fundamental physics. It is based on the tenets of special relativity rather 
than the details of nuclear physics. The speed of sound in a medium, v,, is propor- 
tional to the square root of the ratio of the bulk modulus B (i.e. the ‘strength’) of the 
material to its density p; that is v, oc (B/p)'/*. At supernuclear densities, any material 
that had a high enough bulk modulus to support itself would have a speed of sound 
greater than that of light. As long as relativity theory holds good, this is forbidden. 
This sets a limit in a very fundamental way, because, if the limit were exceeded, there 
could be impossible situations such as an effect preceding its cause (Unit 12). The 
limit set in this way is about 6Mg, but, as mentioned above, the effective limit is 
believed to be less than half this value. 


Having concluded that there is a limit to how massive a neutron star can be, we are 
forced into a corner. Imagine the following thought experiment—of gigantic pro- 
portions. Take two neutron stars, each of mass 2M,, from the thousands in our 
Galaxy, and bring them together so that they begin to coalesce. Their combined 
mass is too great for a neutron star, yet gravity pulls them irresistibly together. 
Continuing collapse will take them well beyond nuclear densities. The strange body 
they would form is called a black hole. Can we predict its properties? 


There are two quite distinct aspects of black hole physics that we will describe in 
this Section. Let us make the distinction from the outset, because the two aspects 
involve very different levels of speculation. First, there is the trapping of light and 
matter (well publicized in science fiction), which we will consider more soberly in 
Section 5.1. This property is predicted by relativity theory, as a direct extrapolation 
of experiments and observations, for any body that is sufficiently compact. Also, it is 
a property that does not depend on the mysteries of super-dense matter. If you could 
fill the solar system with water, then it would trap light for you! Since no radiation 
can escape, such a body is truly black. But having formed such a body (apparently 
unremarkable except for its size), we then encounter the second aspect of black holes. 
Relativity theory also predicts that such a body will start to collapse, with no known 
upper limit for its density. It is here that we face unsolved problems; if we hold 
literally to general relativity theory, the body will inevitably collapse to infinite 
density; and that is difficult to accept. 


General relativity and light 


The relativity theory I have referred to is an extension of the special theory of 
relativity that was discussed in Unit 12. This extension was put forward by Einstein 
in 1915, ten years after he formulated special relativity, and it has come to be known 
as general relativity. You have already seen how the Lorentz transformation of 
special relativity describes the distortion of space and time scales caused by motion. 
Einstein generalized the approach to analyse the way a gravitational field distorts 
space and time scales. Although we cannot go deeply into the ideas involved, some of 
the predictions of general relativity are perfectly straightforward, and in this Section 
we will consider its predictions for the behaviour of light in a gravitational field. 


Gravity affects light in two ways. Both effects contribute to the ‘blackness’ of a 
black hole, and both have been detected and measured, so there is no doubt that 
they exist. These effects are gravitational bending and gravitational redshift, and we 
will examine them in turn*. 


* They are depicted in Course $354, TV 16, so you would find it useful to watch this pro- 
gramme. It is listed in Radio Times under the title Shades of black. 
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Gravitational bending and black holes 


General relativity predicts that rays of light, or of any other electromagnetic radiation, 
are bent when passing through a gravitational field. One way of checking this effect 
is to look for a slight displacement of the image of a star when, due to the Earth’s 
orbital motion, the Sun passes in front of it. The star’s image should return to its 
normal position after a few days, when the Sun has moved well away from the line 
of sight to the star. This effect was first detected in 1919 by a group led by Sir Arthur 
Eddington, and was a triumphant success for general relativity. The bending has 
subsequently been measured much more accurately in several experiments, and the 
same results have been obtained for light and radio waves. 


Measurements have shown that the angle of bending varies inversely as the distance 
of closest approach of the light to the centre of the Sun (Figure 28a), which is in 
agreement with the predictions of general relativity. With the confidence given by 
this success, we may therefore allow the theory to predict what would happen under 
more extreme conditions. Of course, such extrapolation must remain tentative until 
tested by experiment. 


If the Sun were compressed, so that the light grazing the surface passed through a 
region with higher acceleration due to gravity, the bending would be increased. 
Suppose the Sun were compressed sufficiently to become a white dwarf, which may 
well happen some time in the future. Then rays of light can pass closer to the centre 
of the star, and a ray of light grazing the surface, about 6000 km from the centre, 
would be bent through an angle of about three minutes (i.e. 0.05°) (Figure 28b). Now 
throw in another white dwarf so as to make, after coalescence and collapse, a 
neutron star with mass 2M, and with radius about 10 km. Then, if we can rely on our 
extrapolation, the bending near the surface becomes about 70° (Figure 28c). Now 
continue with our thought experiment, and throw in another neutron star. The light 
is then swung right round, and may be trapped in orbit. Here is one reason why such 
an object is called a black hole—it can trap light (Figure 28d). 
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Figure 28 The bending of light by a gravitational field. The angle of bending of a ray of 
light passing by a star (or any other body) is inversely proportional to the distance of the 
ray of light from the centre of the star, and is proportional to the mass of the star. (a) The 
bending near the Sun, as verified by experiment. (b) The bending that would be produced if 
the Sun were compressed to the size of a white dwarf. (c) The bending produced by a neutron 
star with mass 2Ms. (d) The trajectories of light near a black hole. The event horizon is a 
hypothetical surface from within which matter and radiation are unable to escape. 


gravitational bending 
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Gravitational redshift and black holes 


We have seen that gravity deflects light, and this is an indication of one effect of 
gravity on space. But gravity also affects timescales, and the practical effects are 
again quite straightforward. For instance, the frequency of a given source of 
radiation depends on its position in a gravitational field. Thus, identical clocks at the 
base and summit of a mountain should not (and do not) keep the same time. Modern 
atomic clocks are so accurate that this effect can be observed directly. The following 
thought experiment may help to show that there is no mystery in this effect; it fits in 
naturally with other aspects of relativity. 


Suppose I have two tubes, and I suspend identical bodies, such as steel balls, above 
them (Figure 29). I then simultaneously let the balls drop. The ball on the left falls 
through the tube and emerges with kinetic energy Eķin- Applying the principle of 
conservation of energy, this kinetic energy is given by 


Exin = mgH 


where m is the mass of each ball, g is the acceleration due to gravity, and H is the 
height of the tube. For the ball on the right, however, I allow my imagination free 
rein—within the laws of physics. As soon as it drops into the tube, I convert it all 
into radiation; the energy E of the radiation is given by E = mc’. (To achieve this 
would mean starting with a ball made half of steel and half of ‘anti-steel,’ or perhaps 
using an electron and a positron. It does not matter that this is totally impracticable; 
it is possible in principle.) I could beam the radiation down the tube and convert it 
back into a solid ball before it drops out of the end of the tube. But it must drop out 
with the same energy as the ball on the left, otherwise it would be possible to design a 
perpetual motion machine. (For instance, if it dropped out with more energy on the 
right, one could extract some energy from the ball and still bounce it up to the top 
of the left-hand tube.) 


ordinary tube magic tube 


= 


energies measured as balls 
emerge from tubes 


Because the ball on the left gains kinetic energy in falling, we must conclude that the 


radiation also gains energy from the gravitational field as it ‘falls’ down the tube. 
According to the Planck—Einstein formula (E = hf, Unit 13), light that gains energy 
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Figure 29 A thought experiment on 
gravitational redshift. The tube on the 

left is normal. The one on the right is 
imagined to contain a ‘magic’ apparatus to 
change matter to radiation at the top and 
back to matter again near the bottom. 
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increases its frequency. This means it is shifted from the red end towards the blue end 
of the spectrum. Thus falling light is blueshifted. Conversely, light rising from the 
surface of a body is redshifted. As E = hf, the change in frequency Afis related to the 
change in energy AE by the relation 


Af = AE/h 
or alternatively by 
Af AE AE 
FEE 


In the thought experiment, the energy change AE when the radiation rises a distance 
His —mgH, and the total energy E is mc”. So the change in frequency must be given 
by 
Af AE —mgH _ JA 
(se 


(16) 


Now g is equal to GM/r? (Unit 4), where G is the gravitational constant, M is the 
mass of the Earth or star where our thought experiment is being carried out, and r is 
the distance from its centre. Also we may write H, the height of the tubes, as Ar. 
Making these substitutions in equation 16, we get 


Af AE —GMAr 
f E re 


(17) 


Using this result, it is possible to calculate the redshift of a beam of light from the 
surface of a body, like the Earth or a star, all the way out to infinity. The bulk of 
the redshift is attained over a distance of the first few radii of the body, and at a very 
large distance reaches a value 


Af AE GM 
f mae Re? 


(18) 


where R is the radius of the body. The light can lose all its energy (AE = — E) if 
GM/Rc? = 1; that is, if the body has a radius R equal to GM/c?. For an object with 
the same mass as the Earth, this radius is about five millimetres! For a mass equal to 
that of the Sun, however, the radius is about 1.5km, so a neutron star does indeed 
come close to this condition. 


A star that is too massive to be a neutron star will collapse within this limiting 
radius, and thereafter radiation from its surface will lose all its energy before it can 
escape. This is another way of saying that it traps light and becomes invisible*. 


And, of course, if light is trapped, then any form of matter, which always travels 
slower than light, must certainly be trapped. The surface within which matter and 
radiation are trapped is known as the event horizon. For a non-rotating star, it is a 
sphere with radius R given by 

R= ZCM (19) 

E 

The event horizon divides the visible region outside a black hole from the invisible 
region. It corresponds to the surface within which, according to general relativity, 
collapse is inexorable and unending. 


Do black holes exist? 


To answer this question unambiguously in a purely experimental way would involve 
putting a light on the surface of a suitably massive body, stepping back, and watching 
the light go out as the body collapses! Such a pure experiment is, of course, im- 
possible. It also seems hopelessly unlikely that we will ever catch a star in the act of 
collapsing, and record it switching off—even though this must happen if our ideas 
of black holes are correct. However, we can back off a little from such idealizations 
and agree that, if we admit the truth of some basic ideas of relativity, it would be very 
convincing to find a body that has a radius less than 2GM/c?, where M is the mass of 


* Invisible here means not detectable by radiation emitted from the body itself, i.e. non- 
incandescent. This does not preclude detection by being silhouetted against another body, 
or by effects outside the body itself. 


event horizon 
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the body. The mass of such a body might, for instance, be determined by the binary 
star method, even if it is invisible, as was done for Sirius B. But to show that a massive 
invisible body has radius less than a few kilometres could be difficult. Remember 
that if it is as much as 10 km it may be nothing more remarkable than a neutron star. 


There is, however, a ray of hope in the fact that such bodies should not be totally 
invisible. It is only light and matter inside the event horizon that are trapped. In any 
real binary system, there is almost bound to be some matter escaping from the other 
star and being attracted into the black hole (Figure 30). In the process, this matter 
will attain such high speeds before crossing the event horizon that any turbulence in 
its flow may well heat the matter to very high temperatures. Turning this argument 
around, one can say that the only known mechanism for producing very high 
temperatures at the surface of a star (rather than hidden away in the core of a very 
massive star) is turbulent flow of matter towards a very compact massive body— 
either a neutron star or black hole. Calculations show that temperatures of 108 K 
or more can be produced by this mechanism, and probably only by this mechanism. 
At such temperatures, most of the radiation is in the form of X-rays. 


sae 


large companion star 
losing atmosphere into ‘stellar wind’ 


Figure 30 A reconstruction of the supposed structure of Cyg X-1. Cyg X-1 simply means 
the first recorded X-ray source in the constellation Cygnus. The companion to the X-ray 
source is in fact the more massive of the pair. It is assumed that matter from an extended 
atmosphere of the companion star is pulled across to the compact star and is sharply 
accelerated towards it. The stream of gas is heated by violent turbulence, to temperatures 
sufficient to cause the emission of X-rays. 


Unfortunately, to make the delicate distinction between the 10km radius neutron 
star and an event horizon of, perhaps, half that size will be extremely difficult on the 
basis of temperature alone. It will need either a very good model of the system, or a 
fortuitous alignment of the orbital plane of the binary system with the Earth, which 
will lead to eclipses of the compact body that might allow a direct measurement of the 
diameter of the X-ray source. However, if we admit some more relativity theory into 
the argument, we can say that if X-rays are found (indicating temperatures of 
10° K or more) and the mass is greater than 2.5Mg, then there is strong circumstantial 
evidence for a black hole. This is an unsatisfactory argument, in that if we accept 
general relativity in its present form for these extreme conditions in order to give 
significance to a mass that is greater than 2.5M,g, then it follows that black holes can 
exist anyway. The point is rather to look for independent evidence of black holes as a 
test of the validity of general relativity in conditions beyond those under which it has 
already been tested. 


At the time of writing (1981), the most promising candidate for a stellar black hole 
is a compact massive X-ray source, known as Cyg X-/, in a binary pair in the con- 
stellation Cygnus (Figure 30). The masses of the two stars, as estimated from the 
parameters of the orbit, are approximately 20M, and 10Mg. The Jatter is the X-ray 
source, and the X-rays show variations over times as short as 10 ms. Since electro- 
magnetic radiation travels 3000 km in 10 ms, it is argued that the emitting region 
must be less than 3000 km in size. This is much smaller than a star such as the Sun, 
and so tends to confirm the interpretation that the X-ray source is a compact star. 
But it is much too coarse a measurement to make the delicate distinction between a 
neutron star and the event horizon of a black hole of mass 10M,. One cannot talk 
of a ‘theoretical’ radius for a neutron star of mass 10M, because theory says that 
such a star cannot exist, so it is not even clear what distinction one should attempt 
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Cyg X-1 


5.3 


5.4 


to make. This comes from the unsatisfactory nature of an argument that uses general 
relativity to support general relativity. What is needed is an alternative theory, so 
that the characteristics of the Cygnus binary could be used to distinguish between 
them. This point is discussed on the final band of the tape. 


Several other similar objects have been found, but with less convincing estimates 
of the maximum radius. But let us repeat that there is a world of difference between 
saying that they are compatible with being black holes, which is true, and saying 
that they can only be explained as black holes, which is only possible when all 
other interpretations, however outlandish, have been disproved. 


Inside a black hole 


Suppose Cyg X-1 and some of its look-alikes really are black holes, in both senses 
of the term. What happens inside them? Should we accept the prediction of general 
relativity that the density will rise without limit at the centre? Many physicists 
believe that quantum effects will dominate at very small distances in a black hole, as 
they do in the atom. Perhaps we should accept that we cannot describe the centre of a 
black hole until quantum mechanics and general relativity are combined in some 
way to produce a relativistic quantum field theory of gravity. This is where current 
physics can no longer give guidance. Up to now we have been following a path that 
has been laid and trodden by two generations of astrophysicists. Having brought 
you to this point, there are no footprints ahead. 


Now listen to Band 3 of Tape 5, which is about black holes and the far future of the 
Universe—does all matter end its existence in a black hole? 


Summary of Section 5 


1 Neutron Pauli pressure cannot support a star that has a mass greater than about 
2.5M,. There is no known form of matter, and no known force, that could resist 
compression beyond neutron star densities, which are greater than nucleon densities. 
Indeed it appears that, if such matter could exist, it would violate the tenets of 
special relativity by having a speed of sound greater than the speed of light. 


2 Gravity affects electromagnetic radiation in several ways. As predicted by the 
theory of general relativity, radiation is bent by a gravitational field; also its fre- 
quency is observed to be different after having passed through a gravitational field 
from the frequency at emission. Both these effects have been observed experimentally, 
but only for weak fields, which produce very small effects. If general relativity holds 
good for very strong fields, we must conclude that a body of mass M and radius less 
than 2GM/c? can trap light within a sphere, called the event horizon, of radius 
2GM/c?. Once confined within the event horizon, the body cannot resist indefinite 
collapse. This is the definition of a black hole. 


3 Searches for black holes have to contend with two problems. The event horizon 
for a mass Mg has a radius of about 1.5km, which is only slightly smaller than a 
neutron star. In addition, there seems to be little hope of observing the formation 
process of a black hole. Searches have therefore concentrated on binary stars in 
which one member is massive (>2.5M,) and emits X-rays. The X-rays are assumed 
to be generated by matter falling into the massive star. Several such cases have been 
found, but they cannot so far be unambiguously identified as black holes. 


4 (Tape 5, Band 3) As matter is cycled through successive generations of stars, it 
seems that all matter will eventually end up in stars that collapse to black holes. 
This might take hundreds of billions of years, since typical stellar lifetimes are about 
10° years. On a vastly longer timescale, black holes may evaporate by quantum 
tunnelling, and protons may decay. 


Self-assessment questions for Section 5 


SAQ 15 Suppose that you had the power to change the fundamental constants of 
physics, and you wanted to make it possible for a neutron star of mass 10M, to exist. 
What changes would you make? 


SAQ 16 Suppose that the Earth and the planet Pluto were on opposite sides of the 
Sun (Figure 31), and that the Sun were compressed to the size of a white dwarf. Use 
information from Table 1 (Section 4.2) and the data tables in the Appendix to 
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Earth 


Sun asa 
white dwarf 


@ —— 


Pluto 


Figure 31 If the Sun were a white 


dwarf, would Pluto be visible from Earth 
when the planets were on opposite sides 


of the Sun? (SAQ 16). 


determine whether Pluto would be visible from the Earth. (Because Pluto is much 
further from the Sun than is the Earth, you can neglect angle B in Figure 31 com- 
pared with angle a.) 


SAQ 17 (a) If two spheres of water, each of 1 km radius, coalesce and take up a 
spherical shape, what is the radius of the resulting sphere? “(b) If two black holes, 
each with event horizons of 1 km radius coalesce, what is the radius of the new event 
horizon? 


SAQ 18 Calculate the acceleration due to gravity at black hole event horizons of 
radius 1, 10 and 100 km. Is it possible to have a black hole with acceleration due to 
gravity of magnitude 9.8ms ? at the event horizon? 


SAQ 19 What assumptions are made in claiming that Cyg X-1 is evidence for the 
existence of black holes? 


SAQ 20 State which branches of physics are most important in the understanding 
of the following subjects that you have met in this Course. 


(a) The orbit of a planet. 

(b) The Earth’s ionosphere (i.e. plasma at the top of the atmosphere). 
(c) The energy levels of hydrogen. 

(d) The operation, in the future, of a fusion power station. 


(e) Stellar evolution. 


Objectives 


After studying this Unit, you should be able to: 


1 Explain the relation between absolute and apparent luminosity and use the 
relation in simple calculations (SAQs 1 and 3 ; ITQ 1) 


2 Explain the relation between colour (or wavelength) and temperature (SAQs 2 
and 3; ITQ 1) 


3 Plot a Hertzsprung—Russell diagram and identify the regions corresponding to 
main sequence stars and white dwarfs (SAQs 2 and 3; ITQ 1) 


4 Describe how a star is formed from a cloud of gas, and compare the loss of 
gravitational potential energy with the energy radiated from the star (ITQ 2) 


5 Explain the meaning of plasma, and distinguish between electrons, protons, 
neutrons, nuclei, atoms and ions (SAQ 4) 


6 Justify the use of quasi-equilibrium models for stars and discuss how a balance is 
obtained between the various forces and pressures at each point in a star (SAQs 5 
and 8) 


7 Apply the results of kinetic theory of an ideal gas to main sequence stars (SAQs 
8 and 10 ; ITQs 4 and 7) 


8 Explain how fusion reactions occur and how they provide the energy source in 
stars (SAQ 6; ITQs 5 and 6) 


9 Describe the factors governing the rate of fusion reactions in stars (SAQ 7) 


10 Describe the correlations between stellar mass and luminosity, and between 
mass and lifetime (SAQ 9) 


11 Distinguish between the roles of electrons and nuclei in white dwarf stars, and 
distinguish between Pauli pressure and kinetic pressure (SAQ 10) 


12 Relate pressure to the quantum rules governing momentum and the relation 
between momentum and energy (SAQ 10) 


13 Describe how some of the properties of white dwarfs can be deduced from 
observations (SAQs 11, 12 and 14) 


14 Describe how the period and magnetic field of a white dwarf can be measured, 
and estimate the period and magnetic field of a white dwarf, given the period and 
magnetic field of the main sequence star from which it was formed, and assuming 
angular momentum was conserved and flux was frozen in (SAQ 12) 
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15 Explain how the Chandrasekhar limit arises (SAQ 13) 
16 Describe how supernovae and neutron stars are produced (ITQs 10 and 11) 


17 Describe the characteristics of neutron stars and compare them with main 
sequence and white dwarf stars (SAQ 13) 


18 Explain why pulsars were interpreted as rotating neutron stars, rather than 
oscillating white dwarfs or man-made interference (SAQ 14) 


19 Compare the histories of the discoveries of white dwarfs and pulsars, com- 
menting on the nature of the evidence needed to convince scientists of their existence 
(SAQ 14) 


20 Describe how the limits to the mass of a neutron star arise (SAQ 15) 
21 Describe two ways that gravitation affects light (SAQ 16) 


22 Describe the characteristics of black holes, as predicted by the theory of 
relativity (SAQs 17, 18 and 19) 


23 Calculate the size of the event horizon of a black hole of given mass (SAQs 17 
and 18) 


24 Outline a possible way to identify black holes, and comment on the extent to 
which it relies on relativity theory (SAQ 19) 


25 Relate black holes to the very long term future of the Universe 


Appendix: data tables 


Solar system data 


Mean radius of the Earth 6.371 x 10°m 
Mass of the Earth 5.976 x 107+ kg 
Mean radius of the Sun 6.96 x 108 m 
Mass of the Sun 1.989 x 10°° kg 
Period of rotation of the Sun 25.38 days 
Surface temperature of the Sun 5800 K 

Power output of the Sun 3.8 x 107° W 
Moment of inertia of the Sun 5.7 x 10*° kg m 
Solar radiation received at the Earth 1370 W m~? 
Earth-Sun distance 1.5 x 10''m 


General astronomical data 
Distance to nearest star 

White dwarfs 

Neutron stars 

Number of stars in our Galaxy 
Diameter of our Galaxy 
Distance to nearest large galaxy 


z4 x 10! m 

See Table 1, Section 4.2 
See Table 1, Section 4.2 
~ 101 

=~ 50000 light years 
~10’ light years 


Mass and charge of some elementary and atomic particles 


(All numbers quoted are to the nearest integer) 
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mass mass charge mass mass charge 
Particle Mme M, e Particle M, Mn e 
neutrino 0* 0 0 a-particle (i.e. helium nucleus) 7349 4 +2 
electron 1 = =A nucleus with atomic number Z, 
positron 1 = +1 and mass number A = A +Z 
neutron 1839 1 0 
proton } Hepes 1836 1 +1 
deuteron (i.e. deuterium nucleus) 3669 2 +1 


* The mass of a neutrino may not be strictly zero. 
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ITQ answers and comments 


Study comment To make these questions and answers more 
useful for revision purposes, I have noted at the beginning 
of the answer to each ITQ the most important concepts 
required, together with a reference to where those concepts 
are discussed. 


ITQ 1 


Unit 1; Unit 16, 
Section 1.1.2 
Unit 1; Unit 16, 
Section 1.1.2 
TV1; Unit 16, 

Section 1.1.2 


Apparent and absolute luminosity 
Colour—temperature relation 


Hertzsprung—Russell diagram 


To plot the star on the Hertzsprung—Russell diagram, we need to 
know its absolute luminosity and its surface temperature. 


absolute luminosity 


Now apparent luminosity oc 
= , (distance)? 


and so 


= luminosity of star ) z (aen luminosity of = 


absolute luminosity of Sun apparent luminosity of Sun 


( distance of star \? 
ae SSS 
distance of Sun 


1 
= a x (8 x 10°) 
= 0.09 


From Frame 2, a peak wavelength of 800nm corresponds to a 
temperature of 3750 K. When these values of absolute luminosity 
and temperature are plotted on the Hertzsprung—Russell diagram 
in Frame 4, they fall very close to the broken curve. This indicates 
that the star is on the main sequence. 


ITQ 2 
Definitions of energy and power Unit 3 


Total energy output = power x time 
= 107°W x 4.5 x 10° years 
x (no. of seconds in a year) 


The number of seconds in a year is 
3654 x 24 x 60 x 60 = 3.156 x 107 
(It is easy to remember this as approximately n x 107) 


4.5 x 109° x 3.156 x 107J 
1.42 x 10“ 


ll 


So, total energy output 


assuming the mass of the Sun to have been constant. 


1.42 x 10435 
1.989 x 102°kg 


= 7.14 x 10!2Jkg"! 


Energy release per kilogram = 


Gravitational GM, 6.67 x 1071! x 1.989 x 103° = 
energy release | ~ = Ikg 
per kilogram Rs 6.96 x 1 


zx 0.19 x 101? Jkg~! 


Clearly this is much less than the total energy output. (For com- 
parison, the energy released from burning coal or petrol is about 
108 J kg” '.) The simple formula GM,/R only applies to material on 
the surface, which is attracted by the full mass, Mg. Material inside 
will release /ess energy per kilogram. This merely serves to reinforce 
the conclusion that the gravitational energy release falls far short 
of the total energy emitted by the Sun. 
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ITQ 3 
Magnitude of gravitational force = GM, M,/r? Unit 4 
Magnitude of centrifugal force = mv?/r Unit 3 
Mass = volume x density GK 
Volume of a sphere = nr? GK 


Definition of radian Maths in $271 


(a) Fora planet of mass m in a circular orbit with radius r: 


( magnitude of a 


x = (magnitude of the centripetal force) 
gravitational force 


GMm _ mv? 

PSS 
42 

Thus reese 
G 


(b) The star subtends 0.0093 radians, so its diameter is 0.0093r. 
Its mass M is, therefore, 


4n : 4n 0.0093r \> 
M= = Gadus)” x HX = p 


= 42 x 10° ’rp 


(c) The planet travels a distance of 2mr in one year and so its 
speed is 


2nr 2nr A 
t= = Ds 0ra 
lyear 3.16 x 10’s 
(d) From part (b), 
M 
(Sse RS 
42x 10°'r 


Substituting for M using the result of part (a) 
= y? 
P 42x10 7G 


Substituting for v using the result of part (c) 


-eoh Fs 


1.4 x 10?kgm 3 
42x 10°-'G = = 


ITQ 4 


PV = NkT for an ideal gas Unit 11 
Density = mass/volume GK 


For an ideal gas, assumed in Eddington’s model (and most modern 
models) of main sequence stars, 


PV = NkT (eq. 2) 
Thus, N P 
VEE 


Using Eddington’s lower limits, i.e. P = 10'°°Nm~?, T=7 x 
10° K, we get 


N 1045 3 = 
particlesm 
Visca 10S 10° 


= 1.0 x 10°! particles m~? 


Now assuming that the core is composed primarily of ionized 
hydrogen, half the particles will be protons and half will be elec- 
trons. The protons each have mass 1.7 x 10 7’ kg, but the mass 
of the electrons is negligible as far as this question is concerned. 
Thus the density will be 


N = = 
p=m,.x 0.5 x —= 17 x 10-7? x 05 * 1.0 x 107*kgm=? 
V 


=9 x 10?kgm =? 


This is about nine times greater than the density of water. 


ITQ 5 


Ionization energy depends on Z, not A 
E = mc? 


Units 13, 15 
Unit 12 


The ionization energy of the electron in a deuterium atom (i.e. the 
electron’s binding energy in its ground state) is practically the same 
as that for hydrogen, 13.6eV, as the charge on the nucleus is the 
same for hydrogen and deuterium. (There is a tiny difference due 
to the different masses.) The energy released when a deuteron is 
produced in a fusion reaction is obtained by subtracting the mass of 
a deuteron from the sum of the masses of a proton and neutron, 
because a photon has no mass. In terms of the electron mass this is 


(1836 + 1839)m, — 3669m, = +6m, 
Since the mass of an electron is equivalent to about 0.5 MeV, the 
energy released, i.e. the binding energy of the deuteron, is about 
3 MeV. This is more than five orders of magnitude higher than the 
ionization energy. Thus the binding of the nucleus itself is very much 
stronger than the binding of the electron to the nucleus. 
ITQ 6 
(a) Adding all reactions in the first set in a purely formal way gives 
3He + $He + JBe + e` + JLi + p> 
{Be + y + 1.59MeV + JLi + v 
+ 1.37MeV + 23He + y + 17.35 MeV 
Cancelling the particles common to both sides gives 
3He +e + p > $He + 2y + v + 20.31 MeV 
In the same way, adding the second set of reactions gives 
3He + p > $He + 2y + v + e* + 19.28 MeV 
In both cases the nuclei involved are thus 
3He + p > $He 
and about 20 MeV is released. 


(b) The mass of a He nucleus is 7349 m, (from the data tables), 
and this is 


7349 x 9.11 x 10°3'kg = 6.7 x 10 7kg 
The energy release per kilogram is therefore 


20 x 10°eV 


ee a 
- g 


Using the conversion factor 1 eV = 1.60 x 10~!9J, we obtain 
energy per kilogram = 4.8 x 10!4Jkg™! 


Comparing this value with the total energy release of the Sun so 
far (calculated in ITQ 2 to be 7 x 10!2Jkg~'), it is evident that 
nuclear fusion can keep the Sun shining for many times its present 
age. 


ITQ 7 


Meaning of ideal gas Unit 11 


The answer to ITQ 4 was a lower limit of 9 x 10° kgm™? for the 
density of the core of a main sequence star. The average density of 
white dwarfs is higher than this by a factor of about 104, and the 
core density of a white dwarf must be much higher still. Under such 
extreme conditions, matter no longer behaves as an ideal gas. 


ITQ 8 
Energy transferred by a force Unit 3 
Pressure = force/area Unit 11 
Meaning of a derivative Unit 2 


The energy transferred to the gas is the force acting on the piston 
multiplied by the distance moved, because the motion is in the same 
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direction as the force. But force = pressure x area, and so the 
energy transferred is given by 


AE = FAx = PAAx 
Now AAx is the decrease in the volume of the gas in the cylinder, i.e. 
AAx = —AV 


Thus the energy transferred to the gas is 


AE = —PAV 
Th P ae 
35 = SAV 


This expression tells us the average pressure during the volume 
change AV. 


We can express this more precisely using the derivative notation: 


_ limit as (- a dE 


~AV>0\ Aay) 


P SI 
dV 


or, in words, the pressure is equal to minus the gradient of the 
graph of internal energy versus volume. 


ITQ 9 


Inverse square law of gravitation Unit 4 
Energy transferred by a force Unit 3 
Internal energy of a gas Unit 11 


As the mantle moves outwards, the gravitational forces that 
compress it weaken—just as the air pressure is lower at the top of a 
mountain than at its foot. Thus the mantle not only expands out- 
wards, it increases its volume. Initially the expanding mantle is 
heated by the flash of radiation, but once that is over the internal 
energy must decrease as the volume increases, as AE = —PAV 
(ITQ 8). The decrease in the internal energy means that the mantle 
cools. 


ITQ 10 


AE = Amc? Unit 12 


Using information from the data tables in the Appendix, the mass 
difference Am is given by 
Am = m, — (m, + m.) 
= 1839m, — (1836 m. + m.) 
= 2m, 
= 18 x 10° *°kg 


This mass must be produced from the kinetic energies of the proton 
and electron. The minimum amount of kinetic energy needed is 
Exin = Amc? = 1.8 x 1073°kg x (3 x 108ms~!)? 
16-104 
= 1 MeV (as 1.6 x 10° '°J = 1eV) 


ITQ 11 


Rearranging the relation given in the question 
: 7a = ae 1 
c Exin + me? (Exin/mce?) + 1 


l 
= £ x 10° x 1.6 x = i 
9 x 1073! x (3 x 108)? 
= 0.28 


52 
Thus = = 1 — 0.28? = 0.92 
€ 
v = ,/0.92 x 3 ms = 2.9 0 ms 


This is only a few per cent less than the speed of light. 
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SAQ answers and comments 


Study comment To make these questions and answers 
more useful for revision purposes, I have noted at the 
beginning of the answer to each SAQ the most important 
concepts required, together with a reference to where those 
concepts are discussed. 


SAQ 1 
Definition of the tangent of an angle 
Definition of apparent luminosity 


Maths in S271 
Unit 1; Unit 16, 
Section 1.1.2 


The radius of the cluster subtends $° from the Earth, so its length is 


1000 light years x tan (°) = 8.7 light years 
= 8.3 x 10'°m 


(as 1 light year = 9.46 x 1015m) 


The distance of the centre of the cluster is 1000 light years = 
9.46 x 10'8m. 


For nearest member, 


z l > = 
Cesena = x 108 — 8.3 x = — 


= 1.137 x 1073®m~? 


For furthest member, 


1 2 
l/distance)? == 
a ) (sae x 1018 + 83 x ai) sts 
= 1.098 x 10°3®m 2 
j 1.137 
Ratio = —— = 1.036 
1.098 

(The conversion to metres is not necessary; you could use light years 
as a unit if you prefer.) So there is actually a total spread of 3.6% 
in the inverse squares of the distances to the stars. This means that 
the extreme errors caused by using the mean distance to calculate 
the light output from the apparent luminosity will be +1.8%. 


SAQ 2 
R 
white hot 
ee 
cooling 
(a) 
ga 
T 
| white hot 
= 
a7) 
re) 
& 
— n 
2 
cooling 
(b) blue red 
colour 


Figure 32 (a) R versus T for cooling ball. Note that temperature is 
shown decreasing from left to right, as in Frame 4. (b) The corres- 
ponding Hertzsprung—Russell diagram for the cooling ball. 
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Hertzsprung—Russell diagram TV1; Unit 16, 


Section 1.1.2 


Expansion of solid when heated GK 
Colour dependent on temperature Unit 1 
Light output increases (rapidly) with 

temperature Unit 1 


The graph of radius R versus temperature T, and the Hertzsprung— 
Russell diagram, for the cooling iron ball are shown in Figure 32. 
The ball contracts slightly as it cools. We have exaggerated the 
contraction in Figure 32a—the actual contraction is only about 1% 
for a temperature change of 1000 K. On the Hertzsprung—Russell 
diagram, the ball moves towards considerably lower luminosity 
and ‘redder’ colour as it cools. (The total output of radiation varies 
as the fourth power of the temperature, and the light output roughly 
as the third power, but there is no need to remember these details.) 


SAQ 3 

Hertzsprung—Russell diagram TV1; Unit 16, 
Section 1.1.2 
Unit 1; Unit 16, 
Section 1.1.2 
Unit 1; Unit 16, 
Section 1.1.2 


Absolute and apparent luminosity 


Colour—-temperature relation 


Main sequence stars lie around the curve shown on the 
Hertzsprung—Russell diagram in Frame 4. To show that a certain 
star is a main sequence star, we need to plot its absolute luminosity 
versus its temperature (or colour) on such a diagram; if the point 
corresponding to the star falls near to the curve, then we can con- 
clude that the star is likely to be a main sequence star. To determine 
the absolute luminosity, the apparent luminosity and the distance 
of the star must be measured (Frame 1). The temperature of the 
surface of the star can be deduced from measurements of the 
wavelength of the peak in the spectrum from the star (Frame 2), or 
from measurements of the intensity of the absorption lines. 


Alternatively, if the star lies in a cluster, then we can assume that 
the distances of the stars in the cluster are about the same, and we 
can plot a Hertzsprung—Russell diagram of apparent luminosity 
versus temperature for all of the stars in the cluster. This will reveal 
a main sequence curve, and if the star of interest lies on this curve, 
then we would identify it as a main sequence star. 


SAQ 4 
Ionization energy Units 13-15 
Motion of charged particles in magnetic 
field Unit 8 


Atoms are to a certain extent ionized in a Bunsen flame, so that 
electrons and ions are separated. As they are carried upward in the 
flame they are deflected in opposite directions by the magnetic 
field. The flame is therefore widened by the magnetic field, in 
principle if not appreciably in practice. Since the sodium has a lower 
(first) ionization energy than hydrogen, the effect of the magnetic 
field on it will extend to cooler parts of the flame. 


SAQ 5 
Magnitude of centripetal acceleration, 


a= or =— Unit 2 


Magnitude of the gravitational acceleration, 
M E 
g= = Unit 4 
r 


@ = 2n/T Unit 2 


For the Earth, the magnitude of the centripetal acceleration at the 
equator is 

m N 27 2 
a 6.4 x 10°m = | ———__— 6.4 x 10°m 
2E GZ) ES (z x a) = 


= 0.034ms 2 


|centripetal acceleration | 0.034 


0, —— - 3a OS 
| gravitational acceleration | 9.8 


For the Sun, the magnitude of the centripetal acceleration at the 


equator is 


2 2n : 8 
O= 354 dave x 7.0 x_10°m 
.4 days 


2 2 
=( $ ) x 7.0 x 10®ms 2 
25.4 x 24 x 60 x 60 


= 510 = s > 


and the magnitude of the gravitational acceleration is 


GM _ 6.67 x 10°'!Nm’kg~? x 2.0 x 10°°kg 
r? (7.0 x 108m)? 


= 2.7 x 10 ms? 


|centripetal acceleration | 0.0057 
2.7 x 10? 


== z P OS 
| gravitational acceleration | 


Thus the magnitude of the centripetal acceleration on the surface 
of the Sun is rather small and unimportant, being less (relative to the 
gravitational acceleration) than on the surface of the Earth. (Thus 


the Sun is very nearly spherical and, in terms of TV 3, is a flywheel 
that is nowhere near its limiting speed!) 


SAQ 6 
Newtonian approximation for kinetic energy, 
mv? Unit 3 
Conservation of total energy Unit 3 
Meaning of potential energy Unit 3 
z ; 9192 . 
Electrostatic potential energy, E., = Unit 6 
4negr 
Newtonian mechanics valid when v < c Unit 12 


(a) Using the Newtonian approximation, 


Exin = dmv? = 20000eV = 20000 x 1.60 x 10-19J 
= 3.20 x 10715J 


As the mass of a proton is 1.67 x 10° ?’kg, 
1 x (1.67 x 10-27kg) x v? = 3.20 x 10715 J 


v = 2.0 x 10°ms7! 


v 2.0 x 10° 
(b) — = =~—_, = 9.007 

c 3 x 10 
This is much less than 1, so the Newtonian approximation is 
adequate, as assumed in part (a). We could have used the exact 
expression for kinetic energy: 


mc? 


Exin mc? = 20000 eV, 
i-(5) 
e? 


which is a little more difficult to solve, but again gives v = 2.0 x 
106ms 7. 


(c) At the moment of closest approach in a head-on collision 
between particles of equal mass, all kinetic energy is converted to 
potential energy. Thus the potential energy of the system is 2 x 
20000eV = 40000eV. This is positive because the force between 
the particles is repulsive. 


(d) The electrostatic potential energy of two charges of magnitude 
e separated by distance r is 
2 


= 40000 eV 


el 


= 4TEor 
(156 x 10-19)? 
= 
40000 x 1.6 x 10°19 x 4r x 8.85 x 10712 


= 3:6, x 10=14+m 


Thus 


metres 
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SAQ 7 


Strong nuclear force is attractive Unit 15; Unit 16, 
Section 2.2.2 
Coulomb force law Unit 6 
Tunnelling Unit 15; Unit 16, 


Section 2.2.2 
(a) Decrease, because the strong nuclear force is attractive. 
(b) The Coulomb law is 
ES gA }; 
Anegr?\ r 

and this repulsion would have to be made bigger to ‘switch off 
fusion reactions by reducing the probability of close approach. 
This would happen if £ were made sufficiently small. Thus there 


would be no heavy elements produced in a universe with small 
&)—SO presumably no life could exist. 


(c) If the energies of the nuclei were high enough to overcome 
the electrostatic potential barrier anyway, the effect of tunnelling 
would be negligible. Again, if energies were very low, even tunnelling 
would not allow a significant probability of close approach. Of 
course, if one of the particles is neutral (e.g. a neutron), there is no 
electrostatic barrier, so tunnelling becomes irrelevant. 


SAQ 8 
Ideal gas law, PV = NkT Unit 11 
Effect of thermal insulation GK 


(a) The core must be hot enough for the kinetic pressure, which is 
proportional to the temperature for an ideal gas, to be adequate to 
support the star against the (inward) gravitational force. 


(b) The composition affects the transparency of the material to 
radiation. The better the insulating properties of the material 
outside the core the higher the internal temperature. The compo- 
sition of the core has an extra importance in governing the nature 
of the nuclear reactions. 


SAQ 9 


The mass-—lifetime correlation 
Upper limit to white dwarf mass 


Unit 16, Section 2.2.3 
Unit 16, Section 3.5 


(a) The more massive the star, the shorter the lifetime. Hence 
white dwarfs would form first from stars whose mass was about 
1.4 Ms; this is the upper limit—the Chandrasekhar limit—for 
formation of white dwarfs. From Figure 14, such stars have a life- 
time of, very roughly, 2 x 10° years, so that is how long you would 
have to wait! 


(b) As the more massive stars are lost from the main sequence 
earlier than the less massive stars, the average mass of the (re- 
maining) main sequence stars will decrease with time. 


SAQ 10 
Newtonian approximation E,;, = mv? = p?/2m Unit 3 
De Broglie relation p = h/Mag Unit 14 
Relativistic effects occur when v approaches c Unit 12 


(a) Equation 9 relates the pressure of an ideal gas to the mean 
kinetic energy of the particles: 


PV= 3NCExin> 


Using the Newtonian approximation, Ein = p?/2m (where the 
small p represents the magnitude of the momentum), we obtain 


2 
PV = TES 
2m 


and using the de Broglie relation, p = h/Agp, 
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Now there is a difference between <1/A3g> and <Agg>~? (as you 
could easily check by working out the value of each of these averages 
for three values of 4g, e.g. 1 m, 2 m, 3 m). However, for the 


purposes of the remainder of this problem we will ignore this differ- 
ence and write 


Nh? 
PV x ——, 
3mhig 


where ag now just represents a ‘typical’ wavelength for the 
particles. 


(b) From part (a) we can write 


Nh? 1/2 
a eae 
= (r) 


Now (N/V) is the number of electrons per unit volume, which is 
2 x 10°4 electrons per m°. Thus 


z (2 x 1034+ x (6.63 x 10734)? 2 ; 
xX metres 
pe 3 x 9.1 x 1073! x 102° 


= 5.7 x 10°!2m 


= (6.63 x 10734)? : 
2x 91x 10-3! x (5.7 x 10°12)? 


=e ay 


For an ideal gas at 10° K, the mean kinetic energy of the particles is 


<Exiny = 3kT = 1.5 x 1.38 x 107273 x 10°J 
=2x 10!7J 


Thus the typical kinetic energy of electrons in a white dwarf is much 
greater than for an ideal gas at the temperature of the white dwarf. 
We conclude that the electrons in a white dwarf are not an ideal gas. 
In fact as discussed in Section 3.3, the large spacing of the energy 
levels caused by the high density, together with the Pauli principle, 
means that electrons have to occupy states with energies that are 
very much greater than (3/2)AT. 


(d) The typical speed of the electrons can be calculated from the 
kinetic energy: 


dm Ean = 15 X 10S 


£ AE -u = 
So p= [|S ms 
9x 10°33 


1:3¢ 10°ms—* 


& 


This is about half the speed of light. The factor that determines the 
importance of relativistic effects is (1 — (v?/c?)) 1. In this case, 


this factor is 
( (1.3 x a 
1 ek 
(3 x 108)? 


Since this is only 10% greater than the value that applies in the non- 
relativistic limit, I would conclude that relativistic effects are not 
very important. The star would therefore be described as ‘hard’. 


SAQ 1i 
Area of a sphere = 4nr? GK 
Energy loss on cooling = mcAT Unit 11 
Definition of power Unit 3 


(a) The surface area of a sphere is 4nr?, and so the surface area 
would decrease by a factor 


(6.37 x 10° m)? 


Se 
(6.96 x 10ë m)? = 


The absolute luminosity would be reduced by a factor of 8.4 x 1075. 


(b) The power loss of the Sun is about 3.8 x 107° W, so at the 
reduced radius the power loss is 3.8 x 107° x 8.4 x 1075 = 
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3.2 x 10°? W. The energy loss when the ‘collapsed’ Sun cools by 

10* K is 

AE = mcAT = (2.0 x 102°kg) x (4 x 10? Jkg~!K~1) x 10*K 
= 8.0 x 10373 

Thus the time to cool by 10* K is 


8.0 x 10°7J 


(= = 2 5 10'S = 8 x 107 years 
3.2 x 10?7W : 


This is a realistic value for a white dwarfs cooling time. 


SAQ 12 


Conservation of angular 
momentum 

Conservation of magnetic flux in 
highly conducting body 


Unit 4 


Unit 16, Section 3.3.1 


(a) Assuming that no external torques are acting, the angular 
momentum must be conserved. We showed in Section 3.3.1 that 
for uniform spheres this leads to the result 


T; = R;\? 

ENR 
where the subscripts ‘i’ and ‘P denote the initial and final values of 
rotational period and radius, i.e. before and after collapse. 


(eq. 11) 


Thus the period of the contracted star is given by 


7%) = 25.4 = = me 
1 SL a 
‘AR, Tx 108) “7% 


= 1.9 x 10-4 days = 2.7 minutes 


Ti 


(b) We showed in Section 3.3.1 that freezing-in of the field-lines 
in the highly conducting plasma leads to the result 


B; R; z 
FENE 


So the new field strength B; is given by 


R;\? Sy T AOA 
B; = B{—) = 104 x |——_] T 
R; 6 x 10° 


= 14T 


(eq.12) 


As you may have noticed, the values in this question correspond to 
the Sun collapsing to the size of the Earth. 


SAQ 13 


Energy levels in narrow potential 
wells Unit 14 
Relativistic softening Unit 16, Section 3.5 


With such a large mass compressed into such a small volume, the 
energy levels of the electrons would be very widely spaced, and so 
the electrons would have speeds very close to c. The electron gas 
would be too soft to prevent gravitational collapse. 


SAQ 14 (a) Initially, the wavy path of Sirius was thought to 
be due to a companion that was cool (which explained why it could 
not be seen). This was a straightforward, entirely sensible, explana- 
tion. The initial explanation of pulsars was that they were man- 
made interference. This was disproved by showing that the source 
reappeared after one sidereal day. Alternative explanations in 
terms of oscillations of white dwarfs had to be carefully analysed 
and shown to be inadequate. 


(b) First measurements of the temperature of Sirius B were 
dismissed as absurd by some astronomers. Much more recent 
evidence for the temperature, particularly from measurements 
around the peak of the spectrum in the ultra-violet, have led to 
widespread acceptance of the current model of a white dwarf. The 
clinching evidence for the model of a pulsar came from the dis- 
covery of a visible pulsar in the Crab nebula, and from measure- 
ments of the rate at which the period of pulses is increasing. 


(c) A successful model for white dwarfs had to await the develop- 
ment of quantum mechanics by Schrédinger. This came almost a 
century after Bessel’s observations of Sirius. Conversely, neutron 
stars were postulated in the 1930s, soon after the discovery of 
neutrons, and thirty years elapsed before observational evidence 
was presented. 


SAQ 15 


Relativistic effects important when 

v approaches c Unit 12 
The form of the gravitational force GMm/r? Unit 4 
Expression for allowed energies in 


potential well Unit 14 


You could make a more massive neutron star in various ways. A 
large increase in the speed of light would prevent neutrons from 
having relativistic speed, and so increase their Pauli pressure and 
make them harder. A decrease in the gravitational constant would 
decrease the inward force. While you are busy with these, you might 
as well increase the mass of the neutron as well, to slow it down, and 
raise the value of Planck’s constant A, which will raise the values 
of the allowed energies of the neutrons in the potential well (Section 
333): 


SAQ 16 


Trigonometric functions for 
triangle 
Gravitational bending of light 


Maths in S271 
Unit 16, Section 5.1.1 


Pluto would be seen when behind the ‘white dwarf Sun, if the 
bending of light by the white dwarf were greater than «/2 (Figure 
31). 


The angle g is given by 


a radius of white dwarf 6 x 10°m 


2 distance to white dwarf 1.5 x 10°'m 


= 4 x 10 > radians = 2 x 1073 degrees. 


The bending of light by a white dwarf is typically about 3 minutes, 
or 0.05 degrees, which is much larger than «/2. Pluto would there- 
fore never be out of sight. 


SAQ 17 
Volume of a sphere, 42r3/3 GK 
Radius of event horizon, 
R = GM/c? Unit 16, Section 5.1.2 


(a) The total volume of the two spheres of water is 2 x 


(4/3)n(1000)? m3. To make a single sphere with this volume needs a 
radius R given by 


4 4 
T R’ =2 x = (1000) m° 


i.e. R 


1260 m 


(b) A black hole with event horizon of radius 1 km has a mass M 
given by 


Rce? _ 1000 x (3 x 10°)? 


M= = 
G 667 < 10N 
= 1.35 x 10% kg 
Acknowledgements 
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If we take a black hole of twice this mass, M’ = 2.7 x 10°° kg, 
it will have a radius R’ given by 


GM 6-61 K 10 xe ax 102° 
c (3 x 108)? 
= 2000m 


Ra 


Note that the combined black hole has a radius twice that of each 
of the black holes which went into its making. With spheres of 
water, however, the radius only increases from 1000m to 1260m 
when the volume is doubled. In fact we have R oc M for black 
holes, whereas R oc M+’? for spheres of water. 


SAQ 18 
Acceleration due to gravity, 
g = GM/R? Unit 4 
Radius of event horizon, 
R = GM/c? Unit 16, Section 5.1.2 


The acceleration due to gravity at the event horizon of a black 
hole is given by the expression 


where M is the mass of the black hole and R is the radius of the 
event horizon. The radius of the event horizon is given by 


GM 
R = —, and so GM = Rc? 
c 


Substituting this expression for GM into the equation for g above, 
we get 


Thus for the cases given we have 
R= Tkm, g=c?/(1000m) =9 x10 ms? 


R = 10km, g =c?/(10000m) =9 x 10'?ms~? 
R = 100km, g = c?/(100000m) = 9 x 10!!ms~? 


Clearly by making R large enough, g can be brought down to 
9.8ms 7; it just needs R = c?/g = 9.2 x 1015m. The correspond- 
ing mass of the black hole is M = gR?/G = 1.2 x 104° kg, which 
is greater than the mass of the whole Milky Way galaxy! 


SAQ 19 I think the most important assumption is that the 
X-rays are indeed generated by matter falling into a black hole— 
there is no observational proof of this. Next, one must assume that 
the system is a binary (and not a more complicated) system, 
otherwise one cannot deduce the masses of the stars. Thirdly, there 
is the very important reliance on the theoretical result that a com- 
pact star with mass 10 M, cannot be a neutron star. 


SAQ 20 The branches of physics principally involved are: 
(a) Classical mechanics. 

(b) Classical electromagnetism. 

(c) Quantum mechanics and electrostatics. 


(d) Classical electromagnetism, quantum tunnelling, the equiva- 
lence of mass and energy (E = mc’). 


(e) Classical mechanics, electrostatics, quantum theory, rela- 
tivity, nuclear physics. 
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